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The  time  rate  of  change  of  the  potential  energy  is 
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Rayleigh  -  Plesset  Equation 


Vacuum 
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RELATIVE  INTENSIV 


SL  INTENSITY  NORMALIZED  TO  AIR 


WAVELENGTH  (nm) 


SPECTRAL  RADIANCE  (W/nm) 
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cur.  Some  propentci  of  S8SL  inciude  (7)  the  lynchronouj  hollow  cylindrical  P2T  irjn\ducct 
emission  of  iighi  with  each  and  every  acou>iic  cycle,  icin-  was  used  to  drive  the  k’vn.iu.f  m 
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fram  a  RiMi*  WAblii  TW  an  WiMcrMM  gT  ilM  MMNr  of 

ocxamacM  of  a  gl««a  4i,  mwmmi  wyinriily  4Briag  pan  of 
a  babfck  fifaciaML  liw  bia  wUUi  b  flL25  atic  Tfaia  m  h 
«ayf  4«&m4  ai  iha  tfaaa  aT  Um  piavloaa  M  plaa  Hm  ^34 
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laoMatraciad  ftwa  each  af  tba  At  lion  Krta  in  Tigs. 
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Water 


Tcmpcralurt  Dependrner  orSItiele  Bubble  Sunolumlnescence 


•r..i  ^  ’  Delfcf  Lohse.'  wd  William  C.  Mo«s» 

M,hf  „Kl,  ny,A  d,,  Ktmhc/t.  D  J5032  MaHmr,.  Gtn^ny 

U^rmre  Ijvrrmort  ^ioM  Uk^cy.  U^„n«rr.  Califm^  W50 
(RccrivnI  22  Scptcmlier  1997) 

Mic  strtfig  Oependemx  of  the  iMcmky  of  sWIe  botiUt  »cb«  * 

dviunticat  lalculaiMMM  aT  ikm  i  nt  -  ^  .^J[!!""****^*”**^**  approacfi  lo  SflSL  and  futl  ta« 

r»nMcaitaicoiaiHimof!hcbMhMcs»mcfHir  |  $<1031-9007(98 105  VI  I  >9t 


TaCS  mMt«hrr« 

Oi»c  of  (lie  fc3«iir(  of  sfiiple  butihlc 

som,lo,.iirKs.c.Kt  iSBSLl  |I.Z|  «  iIk  w.«i.i.il,  ,.| 
II.C  l.cl.1  t...i«m„  ,u  ,!«,  *.*„  ,<.„^,a«re  e«p<r,i„Kn- 
lally  r.Hind  hy  i|«;  UCLA  fnwp  |2,J|;  tf.  R.  |.  To 
otxa.n  ihcw  >e^ulK.  Baibcf  rl  al  pinmOoi  a,  fo|. 
l-w'  IRcIv  |2.Vi|l:  Wmc.  wa,  eded  l„  ,  ,c.„pc,a. 
lu.c  Ol  2  5  -c  »t.d  c.Hi.plctely  OcgKstii.  An  aif  mes^oic 
ovciliciMi  o(  I  SO  Tihi.  lo  about  209f.  «it  ras 

saioiaiuHi,  iHlj.iMctl  amt  s,H«d.iim«cscciicf  (SI.)  c»- 
ptMimciiiN  unr  pt  iloiu»cd.  Mill  m  2.5  "C.  Ilfcii  iIk  w;»Ici 
«as  iK'aicil  lo  20  X*  Uiitumt  ic^nljuMtug  (lie  pas  c-oiKcn- 
uaiion.  and  ihc  SL  eapciiincni  wa^  repealed.  Finally  the 
same  n,easurcnKrni  was  perforuKd  after  healinp  the  water 
to  32  C.  Ai  all  ilircc  lcm|rcraluie5,  ihe  foreinc  picsMifC 
jhiiiiIhmJc  /*.,  ol  il,c  diiviup  MHiitd  (icld  Wats  adjusted  in 
oidei  lo  Give  masinum,  liphi  iiiicnMiy.  «hi|c  MiaimaiiiiMC 
hoWde  siabilHy  apainM  (ia?iiiciilali(Hi  (u<i/»/f  SL)  Ac- 
i ooliiipioilK-  Mal.il.dl,. Inis'  I*.. 


Strong  huMrte  colbpse  results  tn  a  salisfacliRv  fil  with  lire 
pliysical  values  for  aitd  t'f. 

Because  of  lire  mmplications  in  tire  fits  of  Refs.  12.31 
resulting  data  should  he  read  with  some  care.  This  is 
a^  reflected  in  Fig.  3.  where  we  display  the  data  friHii 
(he«  ti^  icfcrenccs  fm  tlie  c*p.Yn%«m  raiio  (maiinium 
radnw  divided  hy  R„);  if,ey  sliow  large  devialiims 
at  otherwise  tmcltangefl  |tarauicicrs.  11»c  expaiisimi  ratio 
Is  a  quantity  closci)  iclateil  to  tin:  v«,lerKc  of  odhqnc 
and  tlieicloic.  |HeMMM:il.|y.  to  tlie  intensity  of  encigy 
ctmcentmiiun  and  light  eiiiisston  {4|.  ft  is  ihcreftHe 
puuling  that  the  same  light  intensit)  has  been  observed 
in  Refs.  12.31  in  spite  of  tlie  diffcicnt  expansion  ratios 
reported  in  Fig.  3. 

The  cenlial  claim  of  this  pa|rcr  is  ilui  ihc  observed 
wpendcnce  on  w;Hcf  lctn|>craiiite  7  in  Fics  !_3  t;,„ 
be  acciHJntcd  Im  by  ||«:  f  dr,K:odr,Kc  o|  il.e  n.;.iro:.1 


The  acoustic  emissions  from  single-bubble  sonoluminescen< 
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no.).  The  tru.\imum  radiui  i*-l  and  lonoiMfflineiccnce  ^^‘3veleni{(h  tnmt 

■nieniKvi  x )  of  fai  a  nitrnfen  hubbk  and  i  b)  an  arfon  bubble 
are  plotted  o^er  jppfsximaieiy  *40  centeeuitxe  acouMic  eyelet. 

In  boin  easel,  the  acoustic  pmiun  hai  imnally  been  below  the  ,  , 

lijrii*<tnijiti»o  threshold.  The  ptrsture  is  then  rapidly  inertased 
10  a  value  aboie  (he  threshold  I  tame  values  at  in  Fif.  *1. 
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A  new  mechanism  for  sonoluminescence  Lumlntscence  from  SphericaHy  and  Aspherically  Collapsing  Laser  Induced  Bubbles 
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hougl)  ihis  llicor>  lui  been  very  ^u^.'cc^^^ul  m  ihc  vjnc  ol  vjtuui 
jiidonily  cxcacj  (niultiliubbic)  vunuluiiiincseciKC  nccii  hi  The 
iliconc  Oil  i(  ha^  been  mcflicacioui  for  suiioluniiiiCNeencc  clceiin 
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Influence  of  Resonant  rf  Radiation  on  Gas/Liquid  Interface: 

Can  It  Be  a  Quantum  Vacuum  Radiation? 

^  Mirosbv  Collc*-^  *  and  Du  ain  Morse’ 

MD  Om,ionZPM  iJ70  Dnv,  Snilf  C.  CoUla.  CMfomiaVin? 

Mal„u,l,  Dfp,mm,nt.  Ur  t/-.,v<T„n  -./  Cclifomia  a,  S,m,„  B.irhora.  imm  ItmLm.  Calif, ,mia  93106 
(Received  5  Seprember  1997) 

viiT'rhirt^'"". ""  p«x>uc« 

^^nesccncc  h.gh-<2  rewnani  radio  frequency  (tf)  or  microwave  cavity  with  movable  walls  should 

moA-i  I  suggests  that  quantum  vacuum  radiation  might  indeed  be  the  most  feasible 

model  to  cKpiain  both  phenomena,  (S003 1  -9007(98 >05506 *9] 

PACS  nuinben  78.60.  Yj 


The  concept  of  treating  the  gas/wutcr  interface  with 
different  stimuli  to  produce  unusual  effects  has  received 
considerable  attention  in  recent  years  due  to  significant 
interest  in  the  sonoluminesccnce  phenomenon.  Sonolu- 
mincsccncc  (1  -31  is  an  emission  of  light  lesuhing  from 
resonant  ultrasonic  treatment  of  water  containing  bubbles. 
The  generation  of  visible  light  from  24  kHz  ultrasound 
is  an  amazing  amplification  in  frequency  of  1 1  orders 
of  magnitude.  Interestingly,  while  the  dnntiion  nf  iiltra- 
ncumkI  pulses  is  ill  iiiicioscxonds.  the  cmiticd  visilitc  light 
pulses  burst  in  picoseconds  (Ij.  The  mechanisnts  of  the 
sonoluminesccnce  phenomenon  arc  still  largely  unknown. 
Ebcrictn  recently  proposed  a  quantum  vacuutn  radiation 
theory  of  sonoluminesccnce  (41  which  predicts  that  an  os¬ 
cillating  electromagnetic  field  (EMF)  strongly  iiiHucnces 
the  hydrophobic  gas/waicr  interface.  In  (his  model,  the 
two  inicifaecs  witJi  different  Dolartzabiliiies  fwnn-r 


and  surs  ived  freezing,  thawing  cycles  or  boiling  in  u  closed 
container.  It  was  also  realized  that  careful  outgassing  of 
the  water  solutions  resulted  in  the  lack  of  any  measurable 
EMF  effects.  Atomic  hydrogen  seems  to  be  stabilized  in 
a  hydrophobic  hydration  cage  of  argon  or  carbon  Jionide. 
Sonoluminesccnce  phenomena  also  cease  to  exist  in  the 
absence  of  noble  gases  or  carbon  dioxide  (31.  Ouigasscd 
water  has  to  be  sparged  with  gases  containing  either  a 
small  ainotml  i.f  Miihlc  gases  or  caiboii  rlhixhlc  to  poKJucc 
sonoliMiiincscciicc  or  liMI-  dlccts  on  tlic  gas/Jiquid  inter- 
lace.  Our  preliminary  experiments  even  identified  delayed 
rf  emissions  upon  cessation  of  primary  rf  treatment  (un¬ 
published  data).  The  assiKiation  of  two  hydrogen  atoms 
into  a  hydrogen  molecule  radiates  20  cm  wavelength  rfs. 
This  is  the  Irequency  that  radio  telescopes  arc  proeramed 
to  identify.  Others  have  noticed  that  the  presence  o7  water 
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Sensor  Physics:  Signals  and  Noise 

Introduction 

Equilibrium-Thermal  Noise 

Relation  of  fluctuations  to  dissipation 
Total-energy  methods;  frequency  distribution 
Examples 

Shot  and  Nonequilibrium  Noise 
Basic  theory 
Molecular  collisions 
Metals  and  semiconductors 

Sensor  Calibration 

Reciprocity  calibration 
Bessel-null  methods 


Summary 


Given  an  initial  displacement  of 
1  mm,  how  long  will  it  take  for  the 
amplitude  to  decay  to  10-®  mm? 


(A  10-^  mm  amplitude  is 
equivalent  to  an  applied 
acceleration  of  0.5  micro-g’s.) 


Applied  LMbontofy  -  Pann  State 
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^^mMSystems: 


(1-5) 


(1)  times  the  number  of  cycles  required  for  energy  to  decay  by  e  ^  n  times  the  number  of  cycles 
for  the  amplitude  to  decay  by  Alternately,  Q  =  rcN/in(x)  where  A' is  the  number  of  cycles  for 
the  amplitude  to  decay  by  a  factor  of  x, 

(2)  Ratio  of  the  resonance  frequency  to  the  width  of  the  resonance  peak.  The  width  must  be  measured 
as  the  full  width  from  one  half-power  point  to  the  other. 


(3)  For  series -connected  elements:  the  ratio  of  mass  reactance  (or  stiffness  reactance)  at  resonance  to 
the  resistance.  For  parallel -connected  elements:  the  reciprocal  of  that  ratio. 


(4)  In  times  the  energy  stored  in  the  system  divided  by  the  energy  dissipated  per  cycle;  27t  times  the 
resonance  frequency  times  the  stored  energy  di\'ided  by  the  power  dissipated. 


(5)  The  reciprocal  of  2  times  the  damping  factor;  the  reciprocal  of  the  loss  tangent. 

(6)  If  the  damping  is  high  (smalt  Q),  the  resonance  is  isolated  from  other  resonances,  and  there  is  no 
other  mechanism  to  generate  a  changing  phase,  the  Q  can  be  determined  from  the  rate-of-change  of  the 
phase  (in  radians  per  hertz)  at  the  resonance  frequency; 


Q 


fo{df\ 

2  U/J,„ 


MiSystems:  Q 


(7)  From  a  curvc-fu  on  an  HP3562  dynamic  signal  analyzer;  thc.curvc-fit  produces  a  conjugate  set  of  ■ 
poles, ^  +/-  iff,  for  a  resonance  peak.  The  resonance  frequency, and  the  Q  can  be  found  as  follows: 


(1-6) 


/o 

Q 


JIT*  /■  « 


A 

2/, 


where  the  approximations  arc  valid  for  large  Q. 


(8)  Drive  the  system  with  a  square  wave  and  observe  the  ringing  at  the  edge  transitions  of  the  square 
wave.  If  the  peak-to-pcak  amplitude  of  the  first  half-cycle  of  the  ringing  is  a  and  the  peak-to-peak 
amplitude  of  the  second  half-cycle  is  h,  then  the  damping  factor  is: 

_  1  _  \n{a/b) 

22  ^W{alb)'^7t^ 


Note:  The  equivalent  noise  bandwidth  of  a  simple  resonwt  system  is: 
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If  there  is  a  path  along  which  energy  flow  from  a  system  to  its 
environment,  then  energy  from  the  environment  can  flow  back  into  the 
system.  Dissipation  is  a  measure  of  the  energy  that  leaves  the  system  (either 
as  ordered  energy  in  the  case  of  radiation  or  as  disordered  energy  in  the  case 
of  damping);  thermal  fluctuations  are  a  measure  of  the  disordered  energy 
that  enters  the  system  from  the  environment.  In  thermal  equilibrium,  the 
presence  of  dissipation  guarantees  the  presence  of  fluctuations. 


(1-9) 

Total  Energy 

Each  degrec-of'freedom  of  a  system  has  a  **thermar’  energy  of  1/2  kgT  where  kg  is  -  ' 

Boltzmann’s  constant  (1.38  x  lO'^^jouIes/kelvin)  and  7'is  the  absolute  temperature. 

This  thermal  energy  associated  with  each  of  the  components  of  kinetic  energy  {1/2  wv/, 

1/2  ntVy^,  1/2  mv/)^  spring-potential  energy  {1/2  kx^,  rotational  kinetic  {1/2  capacitive 
{1/2  Cf^,  etc. 


A  molecule  in  a  liquid  has 


mv^ 
2  ^ 


A  ball-bearing  in  a  liquid  has 


T 


An  atom  in  a  solid  OR 
A  mass  on  a  spring  has 


-/cx^  =  -icT 

2  2  * 


[The  velocities  and  displacements  indicated  above  are  mean-square  values;  they  represent 
an  average  of  the  actual  fluctuating  velocity  or  displacement.] 


thennl  .pp< 


(I-IO) 


kgT  at  room  temperature:  0.025  eV 

hf  for  visible  light:  1.5  -  3  eV 


acoustic  wave  (100  pPa  in  air): 

acoustic  wave  (100  pPa  in  water): 

chemical  bonds 

covalent: 

ionic: 

hydrogen: 


0.4  eV/cm^ 
30  peV/cm^ 


4eV 

2eV 

0.2eV 


(lev  =  1.6  X  10-'®  joules) 
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magnetite 


(reported  with  penrtssion 
from  Ref.  2) 


magnetotactic  bacterium 

magnetic  dipole  moment,  M  =  1 .3  fA  m^ 
magnetic  flux  density,  B  =  SOptesla 

sensor's  potential  energy  =  MB 
sensor’s  thermal  energy  =  kT 

MB/kT  =  16 


i  Frankel,  Biakemore,  and  Wotfa.  ^WagneWe  in  freshwater  magnetotactic  bacteria.*  Science  203. 1355  (1979) 

2.  Btakamora.  Frankei,  and  KakrtHn.  *Sou1h-sc«king  magnetolactic  bacteria  in  the  Southern  Hemisphefe.-  Nature  2«6.  384  (1980) 


Applied  Research  Labonto/y  -  Pann  State 
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(1-13) 


Distribution  of  Energy 

The  distribution  of  thermal  energy  is  given  by  Nyquist: 

=  ^I^Br^nechanical 

Vl  =  *kiTR,„„„^df 

R  is  resistance;  force  per  velocity  for  mechanical  resistance,  volts  per  ampere  for  electrical 
resistance.  In  general,  the  real  part  of  the  relevant  impedance  is  used  for  R  (which  may  be  a 
function  of  frequency).  df  \%  the  increment  of  bandwidth. 

Since  the  noise  power  is  distributed  over  frequency,  the  noise  is  described  by  a  power 
density  (watts  per  hertz),  or  an  amplitude-squared  per  hertz  (ncwlons-  per  hertz,  volts^  per 
hertz),  or  an  amplitude  per  root  hertz  (pascals  per  root  hertz,  meters  per  second  per  root 
hertz). 

[The  factor  is  about  4x10'^*  at  room  temperature.  These  expressions  are  valid  if  hf « 
kgT  where  h  is  Plank’s  constant.  At  room  temperature,  this  means  that / must  be  much  less 
than  I0»5  Hz.] 


^^^ibnumfriternial  Fluctuations 


spherical  wave  (spherical  source): 

Compute  radial  particle 
velocity  from  Newton’s  Law 
in  fluid: 


p  = 

r 


krj  pc 


Mechanical  radiation 
resistance  (ratio  of  force 
to  velocity): 


Z  .  EA  . 

Ur  l  + 


1  +  (kr)  1  +  (itr) 


Radiation  resistance  for  a  point  r 

source  (A  =  47rr^):  .  ^ 


pc  A{kr^  =  - 


Pressure  fluctuations  associated 
with  “loss”  by  radiation: 


pI  =  =  Akjn^^df 

A  c 


rt  with  a  rigid- walled  box,  side  =  2L, 
sensor  in  the  middle. 

Modes  have  max  pressure  at 
walls  and  at  center: 
cos{lnx/L\  cos(m;^/Z,),  c,os{n7i2/L) 


Wavenumbers  are  then: 
k^  -  Iti/L,  ky  -  rriTi/L,  k.  =  n;i/L 

(spacing  between  k's~  n/L) 
and  k^  —  k^^  +  k^^  + 

Cell  “volume”  in  A:-space  is  (ji/ty  with  one  Z:-point  per  cell. 
Each  mode  gets  kgT {Ml  for  kinetic,  1/2  for  potential)  therefore, 

k-space  density  of  thermal  energy  =  kgT(U7i:)^ 


Divide  by  the  spatial  volume, (2/,/ ,  to  get  the  true  energy  density; 
S  =  471  kj /-’  df/c> 

Another  way  to  write  the  energy  density  is  5  =  p-/ pc^ 
so  the  pressure  fluctuations  associated  with  the  radiation  are  given  by 


Umple  Accelerometer 


(1*2 1; 


R 

-AAAA^ 


(i) 


Vc-v]|^  —j - 


[Impedance  Analogy] 


SYSTEM 

RESPONSE 


0) 

^o> 


2 


^Be'Ut  a  Geophone 


VELOCITY  TRANSDUCER  . 


(1-22 


O.,  = 


mQ 

Cout  =  const.*  {v^-v„) 


=  4k,T-^ 


0)  mQ 


For  the  simple  accelerometer,  (v^  -  is  proportional  to  cd^  below  cDq  and  is  constant 
with  frequency  above  aOg.  The  noise  velocity  (referenced  to  the  case)  is  proportional  to  ar‘. 
Therefore,  the  output  noise  voltage  (1)  is  proportional  to  below  ry^and  to  above  cOg. 


The  equilibrium-thermal  noise  associated 
with  the  electrical  resistance  in  the  sense  coil 
produces  an  output  voltage  noise  (2)  that  is 
independent  of  frequency: 


^out  ^  V  ^  ^^COii 


The  total  noise  is  the  square  root  of  the  sum  of  the  squares  of  the  individual  components. 
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Ultra-Low  Self-Noise  Measurement 

Measurement  of  molecular  agitation  of  32-gram  proof 
mass  in  geophone.  (Approx.  1  nano-g  per  root  Hz) 


^wetf^Optic  Sensors 


(1-24) 


Optical 

fiber 


Noise  associated  with- 

•  Damping  in  spring-mass  system  formed  by 
fiber  and  proof  mass 

•  Optical  scattering  from  thermally  induced 
fluctuations  in  optical  properties  of  fiber 

•  Optical  scattering  from  impurities  and 
density  inhomogeneities  frozen  in  during  fiber 
manufacture 


^ijiiUiBriUm  jPhase  Fluctuations  in  Optical  Fiber 


(1-25) 


L  -  fiber  length 
Aq  =  optical  wavelength 
n  =  index  of  refraction  of  fiber 
a,  =  fiber  radius 

=  optical  mode-field  radius 
6^  =  thermal  penetration  depth  in  fiber 
K  -  thermal  conductivity  in  fiber 


•  Wanser,  et  at..  Optical  Fiber  Sensors  Confencnce. 

Paper  W3.4,  Florence,  Italy,  May  1983. 

•  Wanser.  Electronics  Letters.  2«(1).  53. 1992. 

•  Glenn.  IEEE  J.  Quantum  Electronics  25(6),  1218. 1989. 

r,b<r2  n.t 


The  signal-to-noise  ratio  at  the  output  of  a  linear  circuit 
does  not  depend  on  the  value  of  the  output  load. 


Analysis  of  complicated  circuits  can  often  be  simplified 
by  setting  the  output  load  to  zero  and  calculating  the  ratio  of 
signal  current  to  noise  current. 


The  effective  Q  of  a  system  can  be  changed  by  adding  feedback.  Positive 
feedback  increases  the  Q;  negative  feedback  decreases  the  Q. 

Can  the  noise  of  a  system  be  reduced  by  adding  feedback? 


Not  this  way! 


Shot  Noise  and 
Non-Equilibrium 
Noise 


W^0iestion  #2 


(I!-2) 


•  Measure  the  spectral  density  of 
the  voltage  noise  across  a  10  K 
resistor. 

•  Put  1  mA  DC  current  through 
the  resistor. 

•  By  what  factor  does  the  noise 
voltage  increase?  (Ignore  1/f 
noise.) 


Applied  Reseerch  LMbotmtory~  Penn  State 
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Force 


rate  of  change  in  momentum  of  a  molecule 
initially  traveling  to  the  right  and  hitting  the 
disk  from  behind 


} 


=  (molecular  flux)  (momentum  change  per  collision) 


^0  = 

V  =  2v, 


Looks  like  a  shot-noise  expression! 
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mean-square 

quantity 

average 

fluctuation  in 

=  2 

per 

flux 

flux  density 

carrier 

density 

[j)and\vidtl^ 

[^area^j 


electric-charge  flux  density: 
photon  flux  density: 
momentum  flux  density: 


fn  =  2[^]J(,A//^ 

ll  =  2[hf]l,^flA 

pI  =  2[2/«v]  Pq  A//^ 


Pressure-fluctuation  noise  power  is 
proportional  to  STATIC  PRESSURE 


(II-IO) 


Equilibrium  thermal  fluctuations  in  force  (Nyquist): 


e.g.,  Stokes’ flow  (disk,  radius  a):  P-mech  =  16;;  a 

Pn  =  4  r  I67  a!  A 


Pressure-fluctuation  noise  power  is  almost 
INDEPENDENT  of  static  pressure 


The  shot-noise  foim  requires  that  collisions  be  INDEPENDENT. 
As  long  as  the  mean-free-patfa  is  smaller  than  the  disk  radius,  the 
molecular  collisions  are  highly  DEPENDENT. 

Add  some  kinetic  theory:  Pq  ~  n  kg  T 

7]  nmv(rrifip)/3 

(PniV  /  (Pn2^J  3  7c/ 8  (radius)  /  (mean-free-path) 


smmciors 


(11-12) 


FERMI 

LEVEL 


Carriers  (electrons)  are  highly  correlated 
Noise  is  independent  of  flow  volume  (current) 

=  (4t,r/«)A/ 


FERMI 

'  level 


Carriers  (holes  or  electrons)  are  independent 
Noise  is  dependent  on  flow  volume  (current) 
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Energy  (eV) 


Applied  Voltage  (n-i9) 


For  qAV  «  kj  2^  - = - \df  =  4k,T^ 

R  {qAV/k„Tj  “  R 


For  qAV  «  kgl 
For  qAV  »  k^T 


Example  from  measurements  hi 
search  of  fractional  quantum 


Samtnadayar.  Glaitli,  /in,  and  Etienne, 
Ph>-s  Rct'.  Lett.  79,  2526  ( 1997) 


2^(1)#  =  Iqldf 


Shot  Noise 
for  q/3 


0  100  200  300  400  500  600 


J 

h 

mmoise 

For  example,  if  the  current  is  produced  by  random  events  and  each  event  has  an 
exponential  decay  (with  time  constant,  tg): 

oo 


^(0  =  ;  lit)  =  e-"‘>u(t) 

/  =  ! 


then  the  transfer  function  in  the  frequency  domain,  H(jo)X  gives  the  spectral  shape 
of  the  noise  power: 


If  there  are  many  processes  that  can  be  triggered  by  the  random  impulses  and 
each  process  is  exponential  with  its  own  unique  time  constant,  then  the  spectral 
distribution  of  the  noise  power  can  depart  significantly  from  either  white  noise  or  a 
I/f  power  distribution.  This  may  be  the  way  1/f  noise  distributions  are  produced. 
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fj[-^^6is^e:'^ljBrief  land  Inadequate  Introduction  (ii-23) 


MANY  physical  processes  produce  fluctuations  with  a  power  spectrum 
that  goes  as  1/f. 

The  noise  power  in  excess  of  the  equilibrium-thermal  fluctuations  is 
associated  with  power  input  to  the  system  that  drives  the  system  away 
from  equilibrium. 

Observed  1/f  noise  can  extend  over  many  decades  in  frequency.  If  the 
multiple-exponential-process  model  is  correct,  then  there  must  a 
correspondingly  large  spread  in  process  time  constants. 

The  integral  over  all  frequency  of  a  1/f  power  distribution  is  infinite  so 
there  must  actually  be  a  lower  limit  to  the  Iff  behavior.  This  means  that 
there  are  at  least  two  free  parameters:  the  total  fluctuation  power  and  the 
lower  frequency  limit. 
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Summer  School  - 1 998 


(iii-i) 


Sensor 

Calibration 


(111-2) 


-a 


e  =  A  i  ^  B  V 
f  ~  C  z  4-  D  V 

If  5  =  C  then  the  device  is  reciprocal: 


Reciprocal  transducers: 

electrodynamic  (moving-coil) 

piezoelectric 

capacitive  (small-signal) 


Nonreciprocal  transducers: 
piezoresistive 
electron-tunneling 
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'^j^^^rocity 


(111-3) 


If  B  =  -C  then  interchange  the  roles  of  the 
variables  on  one  side  of  the  device. 

/ 

- ► 

+  0 - 

e 

-o - 


■0“h 

e  =  A  /  +  B/ 

V  =  C  /  +  D/ 

■o- 


A  =  (AD-CB)/D  B  =  B/D 

C  =  -C/D  D  -  1/D 


Therefore  B  =  C 


If  B  does  not  equal  either  C  or  -C.  then  the  device  is 
nonreciprocal. 


If  a  device  is  reciprocal,  then  the  ratio  of  the  output  potential  to  the 
inputyZow  is  the  same  regardless  of  which  port  is  taken  to  be  the  input. 

For  example: 


(^1  ^2)  h  ^2  ^2 

^2  =  ^2  h  ^2  h 


~  “  ^  (A  transfer  impedance  in  general.) 

h  /.=0 


Reciprocity  does  NOT  depend  on  the  device  being  lossless. 


Reciprocity  does  A^OTmean  that 
voltage  ratios  are  identical. 


10  V 

dQ.v 
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(111-5) 


Given  a  transducer  (not 
necessarily  reciprocal): 

What  is  its  response? 

Receiving  response: 


V 


a 


Transmitting  response:  P  ~ 


It  is  often  inconvenient  and  inaccurate  to  directly  measure 
forces  and  velocities. 


lalibration 


(III-6) 


Use  three  transducers:  a  source  (S),  a  receiver  (R),  and  a  reciprocal 
transducer  (T).  Altogether  there  are  four  unknowns:  a^,  aj,  pj,  and  p^. 

Connect  the  transducers  through  a  known  transfer  impedance  and 
make  the  following  measurements: 


^S1 


xfer 


R 


^xftr 


•rs 


■^xfer 


R 


o 

^Rl 

■o 

[SRI 

o 

■o 

[ST]  =  ej^2  /  ^2 

o 

[TR]  = 

This  gives  three  equations  for  the  four  unknowns.  Reciprocity  provides 
the  fourth  equation  and  allows  solving  for  all  four  responses. 
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a)  Measure  currents  by  measuring  voltage  across  a  resistor.  Resistance  and 

voltage  ratios  are  easier  to  measure  accurately  than  absolute  voltages. 

b)  Set  isi  ~  ^2  and  adjust  so  that  This  produces  the  same  field 

at  the  receiver  location  for  each  measurement  and  also  permits  the  use  of 
source  and  receiver  that  are  not  linear. 

c)  If  two  reciprocal  transducers  are  available,  measure  the  reciprocity 

explicitly  to  check  the  transducers  and  apparatus. 

d)  If  two  “identical”  reciprocal  transducers  are  used,  then  only  two 

measurements  are  required. 

e)  In  some  circumstances,  only  one  transducer  is  required:  (1)  excite  a 

lightly  damped  system  and  measure  response  during  decay,  (2)  transmit  a 
pulse  and  measure  a  reflection. 


Traveling-Wave  Tube 


^^^sel-NUll^  Calibration 


LASER 


EMMM^alibration 


FIBER 


Component  that  reflects  from  moving  mirror: 

[1]  =  Acos{Q)Qt +  2kQd) 

Component  that  reflects  from  cleaved  end  of  fiber: 

[2]  =  B  cos  {co^t) 

Photodetector  output  (square-law  detector): 

PD  = 

PD  cos(2 

Sinusoidal  motion  of  mirror; 

d  -  Jq  +  sin(ty,f) 


Adjust  the  drive  level  to  null  the  output  of  the  photodetector  at 
the  drive  frequency.  These  nulls  correspond  to  the  zeros,  Zy,  of 
the  Bessel  function,  Jj.  Consequently,  the  displacement 
amplitude,  df,  is  only  a  function  of  the  laser  wavelength, 


d 


1 


_fL_ 
2  k, 


An 


(where  Z;  =  3.83171,7.01559,  10.17347,  13.32369,...) 


(III-21) 


Sensor  System 


desired 

signal 


■  ■.-a'iP'OST- 

PROGESSING 


Acoustics  Summer  School  - 1998 

, 

Sensor 

Electronics 

Supplement 

Often,  but  certainly  not  always,  the  noise  floor  of  a  sensor  system  is  set  by 
the  first  stage  of  electronics  connected  to  the  transducer.  In  the  design  of  high- 
performance  sensors,  it  is-important  to  understand  the  interaction  between  the 
sensor  and  the  electronics. 
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OBSERVATIONS- 

Noise  at  the  output  of  a  preamplifier  depends  on 

(1)  frequency 

(2)  impedance  of  the  sensor 

(3)  gain  of  the  preamplifier 

SIIVIPLIFIED  MOOFI  • 

Eliminate  (3)  by  referring  noise  to  preamplifier  Input 

Simplify  dependence  on  (2)  by  using  equivalent  voltage 
and  current  sources 

Applied  Research  Laboratory  -  Penn  State 

foni  1 1 


In  order  to  understand  this  interaction,  the  characteristics  of  preamplifier 
must  be  examined.  If  a  measurement  of  output  noise  from  a  preamplifier  is 
made  with  various  sensor  elements  connected  to  the  input,  the  output  noise  is 
found  to  depend  strongly  on  three  factors:  frequency,  the  impedance  of  the 
sensor  connected  to  the  input,  and  the  gain  of  the  amplifier. 


In  many  cases,  the  output  noise  is  directly  proportional  to  amplifier  gain 
(suggesting  that  the  dominant  noise-producing  mechanisms  in  the  amplifier  are 
in  the  input  stage)  so,  by  referring  the  noise  to  the  amplifier  input,  the  effects 
of  amplifier  gain  can  be  eliminated. 


Also  (in  most  cases)  the  noise  can  be  divided  into  two  components;  one  that 
is  independent  of  the  impedance  of  the  attached  sensor  and  one  that  depends 
linearly  on  the  magnitude  of  that  attached  impedance.  Consequently,  the 
amplifier  can  be  well  represented  by  an  equivalent  voltage  noise  and  an 
equivalent  current  noise  (the  current  noise  being  the  component  that  produces 
the  noise  proportional  to  the  attached  impedance). 
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If  the  noise  of  an  amplifier  is  carefully  measured  with  many  different 
resistors  connected  to  the  input  (to  simulate  a  wide  range  of  sensor 
impedance),  the  results  appear  as  shown  above.  Below  some  value  of 
resistance,  the  noise  is  constant.  This  value  is  selected  for  the  value  of  the 
equivalent  voltage  noise  component.  Above  some  higher  value  of  resistance, 
the  noise  increases  linearly  with  resistance  and  the  equivalent  current  noise 
value  is  obtained  from  the  slope  of  that  line. 


The  resistor  produces  noise  of  its  own  (from  equilibrium  thermal 
fluctuations  in  the  material  —  Johnson  noise)  and  for  a  good  low-noise 
amplifier  (A2),  that  noise  ciin  be  measured  directly  over  some  range  of 
resistance.  A  poor  amplifier  (Al),  on  the  other  hand,  shows  no  such  region. 
The  closest  a  poor  amplifier  comes  to  reaching  the  resistor’s  Johnson  noise  is 
for  the  value  of  resistance  equal  to  the  equivalent  noise  voltage,  e^,  divided  by 
the  equivalent  noise  current,  The  region  over  which  a  good  amplifier  can 
measure  resistor  Johnson  noise  extends  symmetrically  about  the  value  e^/i^. 
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To  analyze  the  noise  performance  of  a  sensor-amplifier  combination,  the 
noisy  amplifier  can  be  replaced  by  a  noiseless  amplifier  with  a  noise-voltage 
source  and  a  noise  current  source  attached  to  the  input.  (The  noiseless 
^plifier  still  has  the  same  input  and  output  impedances  as  the  real  amplifier, 
though.)  For  most  purposes,  it  is  acceptable  to  treat  these  two  sources  as 
completely  uncbrrelated  even  though  there  may  be  some  correlation  in  reality. 
In  the  case  of  complete  correlation,  the  maximum  error  that  can  be  introduced 
is  a  factor  of  the  square  root  of  two  in  amplitude. 

In  general,  there  is  some  frequency  dependence  associated  with  these  noise 
sources.  Both  the  voltage  noise  and  the  current  noise  can  have  a  low- 
frequency  dependence  of  1/f  in  power.  In  addition,  the  current  noise  often 
increases  with  frequency  above  some  frequency.  Components  and  amplifiers 
intended  for  low-noise  applications  will  normally  have  voltage-  and  current- 
noise  spectra  provided  by  the  manufacturer. 
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By  selecting  and  i^,  a  two-dimensional  space  can  be  constructed  on 
which  various  sensor  materials  and  elements  can  be  compared  to  various  types 
of  amplifiers.  A  number  of  off-the-shelf  amplifiers  and  discrete  transistors 
suitable  for  sensor  preamplification  are  shown  on  this  diagram  (with  the  noise 
evaluated  at  1  kHz). 

JFET  devices  are  especially  suited  to  low-noise  applications  in  which  the 
current-noise  component  must  be  minimized;  BJT  devices  are  particularly 
suited  when  it  is  desireable  to  minimize  the  voltage-noise  component.  (The 
line  coimecting  three  of  the  symbols  illustrates  the  performance  of  a  single 
transistor  under  different  conditions  of  collector  current.)  MOS  devices  are 
useful  at  ultrasonic  fi-equencies  and  beyond  but  are  plagued  with  high  1/f-noise 
powers  that  can  dominate  performance  at  low  fi-equency.  In  addition,  several 
micropower  (ppwr)  devices  are  included  since  overall  power  consumption  is 
often  a  critical  specification  in  a  sensor  system. 

It  is  readily  apparent  that  it  is  difficult  to  achieve  much  better  than  1 
nanovolt  per  root  hertz  or  0.5  femtoamp  per  root  hertz  regardless  of  device 
type.  Of  particular  note  is  the  fact  that  the  space  is  not  covered  uniformly  by 
devices.  In  general,  cost  and  power  consumption  increase  in  the  direction  of 
decreasing  voltage  noise. 
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When  connected  to  a  sensor,  the  internally  generated  noise  can  be  described 
by  three  components:  (1)  the  amplifier  voltage  noise,  (2)  the  amplifier  noise 
current,  flowing  through  the  sensor  producing  a  noise  voltage  equal  to  that 
noise  current  times  the  magnitude  of  the  sensor  impedance,  and  (3)  a  thermal- 
equilibrium  (Johnson)  noise  component  associated  with  the  real  part  of  the 
sensor  impedance.  The  plot  shown  above  is  representative  of  a  piezoceramic 
sensor  for  which  the  sensor  impedance  is  primarily  capacitive  and  for  which 
the  real  part  of  that  impedance  is  dominated  by  dielectric  loss.  Not  shown  in 
this  diagram  are  external  resistors  for  bias  and  gain  as  would  be  used  with  an 
op  amp.  These  resistors  produce  Johnson  noise  and  also  interact  with  the 
amplifier  noise-current  source  but  the  accounting  is  straightforward.  Each 
noise  source  is  considered  in  isolation  and  the  results  are  root-mean-square 
summed. 


For  any  sensor,  it  is  crucial  to  distinguish  between  the  magnitude  of  the 
impedance  and  the  real  part  of  that  impedance  since  the  magnitude  determines 
the  effects  of  amplifier  current  noise  (and  so  is  a  function  of  amplifier 
selection)  while  the  real  part  produces  a  noise  component  unrelated  to  the 
amplifier. 
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The  irreducible  component  is  equilibrium-thermal  noise  internal  to 
the  sensor  and  depends  on  the  resistive  (real)  part  of  the 
sensor  impedance,  . 

M  is  the  sensor  response:  volts  per  pascal  in  this  example. 

The  amplifier  contributes  two  terms,  one  of  which  depends  on  the 
sensor  impedance. 


ApplM  Res«*rch  Laboratory  -  Pann  State 


A  useful  form  in  which  to  examine  the  various  noise  corhponents  is  in  terms 
of  noise-equivalent  pressure.  If  the  voltage  response  (volts  per  pascal,  for 
example)  of  a  pressure  sensor  is  M,  then  voltage-noise  terms  can  be  referred  to 
equivalent  pressure  by  dividing  by  M.  To  properly  express  the  incoherent 
addition  of  noise  terms,  the  expression  above  is  written  in  mean-square 
pressures  and  voltages.  There  is  an  ambient  noise  term  a  term 
connected  with  the  internal  loss  in  the  sensor  (hence,  “irreducible”),  and  two 
terms  associated  with  the  amplifier. 


If  the  noise  floor  is  dominated  by  the  amplifier  contribution,  there  are  two 
terms  to  consider.  In  a  particular  situation,  one  of  those  terms  may  be 
considerably  larger  than  the  other.  However,  if  it  is  not  known  which  (if 
either)  term  dominates,  then  a  compromise  term  can  be  introduced  as  the 
geometric  mean  of  the  two  amplifier  terms.  Clearly,  this  is  not  a  good  strategy 
if  one  of  the  terms  is  much  larger  than  the  other.  In  the  special  case  in  which 
the  two  amplifier  terms  are  of  the  same  order,  the  geometric-mean  term  is 
useful. 
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In  designing  a  high-performance  sensor,  the  desired  noise  floor,  often 
determined  by  the  ambient  background,  should  be  identified.  When  this  is 
coupled  with  the  design  bandwidth  of  the  sensor ,  the  noise  components  can  be 
individually  assessed  against  this  noise-floor  specification.  Ignoring  for  the 
time  being  the  term  resulting  from  dielectric  loss  (term  2),  the  design  for  a 
prezoceramic  hydrophone  with  a  self  noise  below  quiet  ocean  ambient  in  the 
band  10  to  1000  Hz  might  appear  as  shown  above.  At  the  lower  band  edge,  the 
desired  backgroimd  sets  an  upper  limit  on  the  allowable  amplifier  current 
noise,  while  at  the  upper  band  edge,  the  background  sets  an  upper  limit  on  the 
allowable  amplifier  voltage  noise.  The  respective  amplifier  components  can, 
of  course,  be  lower  than  these  limits.  This  limiting  configuration  can  be 
specified  in  terms  of  the  intersection  between  components  3  and  4  given  by  p. 
and/o. 


76 


If  a  capacitive  sensor  were  being  designed  for  an  ultrasonic  band  (100  kHz 
to  10  MHz),  then  the  maximum  permissible  levels  of  voltage  arid  current  noise 
change  as  shovm  above.  The  principle  is  still  tiie  same.  Identify  the  desired 
noise  floor  specification.  Overlay  the  voltage  and  current  components  of 
amplifier  noise  adjusting  e^  and  i^  until  the  hypothetical  amplifier  is  just 
adequate.  A  real,  acceptable  amplifier  would  then  be  one  for  which  the  voltage 
and  current  noise  components  would  be  equal  to  or  less  than  these  values. 
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For  a  sensor  having  an  impedance  that  is  primarily  resistive  (a  geophone  or 
a  piezoresistive  pressure  sensor,  for  example),  the  amplifier  terms  (to  a  first 
approximation)  are  constant  with  fi-equency.  The  same  analysis  applies:  both 

amplifier  terms  must  be  at  or  below  the  desired  noise  floor  over  the  design 
bandwidth. 
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preamplifier  Performance  Measures 


Noise 

Noise  Power 

“resistance" 

Coefficient 

R  =  ^ 

II 

i 

A 

h 

4kJ 

JFET 

1M- 100M 

0.0001  -0.001 

BJT 

1K-1M 

0.  01  -0.  1 

Applied  Research  Laboratory  •  Pertn  State 


There  are  two  useful  measures  of  amplifier  noise  performance.  A 
commonly  used  parameter  is  the  noise  “resistance”  given  by  the  ratio  of  noise- 
voltage  spectral  density  to  noise-cuitent  spectral  density.  This  is  often  cited  as 
the  value  to  which  the  resistive  part  iof  the  sensor  impedance  should  be 
matched  for  optimum  noise  performance.  As  discussed  below,  there  is  some 
truth  to  this  statement  but  the  merit  of  an  amplifier  is  hot  determined  by  how 
close  its  noise  resistance  is  to  the  source  resistance. 


Another  important  (but  little  used)  measure  is  the  noise  power  coefficient 
defined  as  the  ratio  of  the  product  of  the  voltage-  and  current-noise  densities  to 
the  thermal-noise  power,  4kgT.  This  is  a  better  measure  of  the  performance  of 
an  amplifier  than  either  the  noise  voltage  or  noise'  current  alone.  The  smaller 
this  ratio,  the  better  the  amplifier  is.  However,  the  noise  power  coefficient  is 
not  sufficient  to  determine  the  overall  performance  for  a  particular  sensor. 

Both  of  the  above  measures  are  necessary.  (The  analysis  of  an 
amplifier/sensor  combination  can  be  performed  perfectly  well  using  both  the 
current-  and  voltage-noise  spectral  densities  without  reference  to  these  other 
measures  but  there  are  some  conceptual  advantages  to  the  measures  given 
above.) 
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In  general,  the  noise  floor  should  be  controlled  by  elements  close  to  the 
system  input.  For  many  designs,  proper  selection  of  the  amplifier  can  result  in 
the  amplifier’s  noise  floor  being  below  the  floor  set  by  theisensor  itself  The 
first  condition  above  applies  to  the  case  of  a  purely  resistive  sensor.  The 
smaller  a  is,  the  larger  the  range  of  source  resistance  is  that  can  be 
accomodated  without  masking  by  a  particular  amplifier. 


If  the  sensor  impedance  is  not  purely  resistive,  then  the  second  condition 
above  specifies  the  ranges  of  parameters  for  which  the  amplifier  will  not  mask 
the  sensor  noise.  This  case  is  more  complicated  because  the  impedance  is,  in 
general,  a  function  of  frequency  so  the  conditions  may  be  satisfied  over  some 
regions  of  the  spectrum  and  not  satisfied  over  other  regions.  In  a  number  of 
important  cases,  the  ratio  of  the  resistive  part  of  the  impedance  to  the 
magnitude  of  the  impedance  (the  “loss  tangent”)  is  approximately  constant 
over  a  range  of  frequency,  which  simplifies  the  determination  somewhat. 


In  either  case,  if  a  (or  a’)  is  less  than  one,  there  is  a  region  over  which  the 

sensor  noise  (the  component  from  the  real  part  of  the  impedance)  MAY  set  the 
noise  floor. 
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ypreamplifier/Sensor  Interaction  the  irreducible  term 

Real  part  of  the  sensor  impedance: 

l-^sl 

1.  Purely  resistive  sensor:  ^  = 

Re[Z,] 

1 

2.  Primarily  reactive  sensor  (e.g.  capacitive): 

\^S\  ^ 

RelZ,]  ■  ^ 

3.  Typical  sensor: 

•  electrical  loss  (loss  tangent,  S) 

•  electrical  loss  (resistance.  R^) 

•  mechanical  loss  {Q„) 

•  radiation  resistance 

Applied  fiesetrch  Lebonlory  •  Penn  State 

•*»l  1 

The  ratio  of  impedance  magnitude  to  real  part  is  one  for  a  purely  resistive 
sensor  and  is  equal  to  the  reciprocal  of  the  loss  tangent  for  a  reactive  sensor. 
Since  the  electrical  impedance  of  the  sensor  includes  the  electrical  equivalent 
of  the  sensor’s  mechanical  elemerits  (by  means  of  the  transduction 
mechanism),  the  losses  and,  therefore,  the  noise  may  be  caused  by  either 
electrical  losses  or  mechanical  losses. 
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Applied  Research  Laboratory  •  Penn  State 


A  number  of  various  types  of  amplifiers  (either  op-amps  or  discrete 
transistors)  are  plotted  above  on  the  coordinates  of  noise  power  coefficient  and 
noise  resistance.  This;chart  is  useful  for  selecting  an  amplifier  type  once  the  > 

sensor  impedance,  the  operating  bandwidth,  and  the  desired  noise  floor  are 
fcnovim. 
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ij^W&Exiample 


Geophone  array  formicroseismic  measurements 


Performance  with  respect  to  seismic 
acceleration  background,  g: 


One  Geophone 


^©quiv 

[nV  prH] 


One  Hundred  Geophones 
in  Series 


Applied  Ras^Mfch  Lmboratory-  Peon  State 


•oKjiJX 


This  is  an  illustration  of  the  design  of  a  geophone  array  for  sensing 
microseismic  disturbances.  The  frequency  range  of  interest  is  0.2  to  2  Hz.  The 
curve  labelled  gM  is  the  specification  for  self  noise  for  this  sensor.  (This  curve 
is  produced  by  multiplying  the  ndise-floor  in  seismic  acceleration  by  the 
trahsfer  fimction  of  the  geophone  or  geophone  array.  This  produces  an 
equivalent  voltage-noise  noise  floor.)  For  a  single  geophone,  both  the 
amplifier  voltage  noise  and  the  irreducible  noise  associated  with  the  electrical 
resistance  of  the  geophone  coil  are  well  above  the  desired  noise  floor.  By 
connecting  100  geophones  in  series,  however,  the  desired  noise  floor  is 
achieved.  This  may  not  be  an  elegant  solution  but  the  result  is  cheaper  and 
more  rugged  than  the  closest  commercial  alternative. 
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Amplifier  Noise  Specifications 


Equivalent  Voltage  and  Current  Noise  Values  for  Common  Amplifiers 

T~  '  [  1  —  -7 -  . .  ,,  - _ 


j 

1 

> 

1  iA 

1 . 

eA/iA 

i 

eAiAMkT 

nVprH 

fAprH 

ohms 

! 

BJT 

LT1007 

2.5 

400 

6250 

0.0625 

- 

Lt1012 

22 

70 

314286 

0.0963 

LT1024 

14 

20 

700000 

0.0175 

-  - -  . .. 

LT1028 

0.85 

1000 

850 

0.0531 

- - - - 

OP05 

10 

140 

71429 

0.0875 

..  - - 

OP27 

3 

1000 

3000 

0.1875 

OPA77 

7.5 

220 

34091 

0.1031 

— 

CLC425 

1.05 

1600 

656 

0.1050 

{1mA] 

MAT02 

0.85 

750 

1133 

0.0398 

- - 

(luA) 

MAT02 

20 

25 

800000 

0.0313 

— 

(1mA) 

LM394 

1 

500 

2000 

0.0313 

(luA) 

LM394 

20 

25 

800000 

0.0313 

- - 

“ 

JFET 

LT1022 

14 

1.8 

7777778 

0.0016 

LT1055 

15 

2 

7500000 

0.0019 

LF353 

16 

10 

1600000 

0.0100 

”  - 

LF411 

25 

10 

2500000 

0.0156 

— -  - - - 

_ _ 

AD743 

3.2 

6.9 

463768 

0.0014 

_ 

AD744 

18 

10 

1800000 

0.0113 

- - - 

AD549 

35 

0.2 

-175000000 

0.0004 

— 

_  _  _  . 

OPA111 

8 

i  0.6 

13333333 

0.0003 

- - - - - 

U401 

2 

0.5 

4000000  ‘ 

0.0001 

2N4338 

6 

0.8 

7500000 

0.0003 

2N6485 

7 

1 

7000000 

0.0004 

uPwr 

LM4250 

50 

80 

625000 

0.2500 

OP20 

60 

80 

750000 

0.3000 

OP21 

20 

200 

100000 

0.2500 

OP22 

90 

180 

500000 

1.0125 

OP90 

60 

700 

85714 

2.6250 

OP191 

35 

800 

43750 

1.7500 

OP193 

65 

50 

1300000 

0.2031 

OP196 

26 

190 

136842 

6.3088 

OPA1013 

25 

120 

208333 

0.1875 

MOS 

CA3440 

110 

1.8 

61111111 

0.0124 

CA3140 

40 

1.8 

22222222 

0.0045 

CA3160 

72 

0.8 

90000000 

0.0036 

(chopper) 

LTC1052 

30 

0.6 

50000000 

0.0011 

Notes 

1.  eA  and  i/ 

are  the  equivalent  voltage  and  current  noise  components  respectivelv 

r 

iVprh  is  nanovolts  per  root  hertz;  fAprh  is  femtoamps  per  root  hertz  T 

1 

All  values  are  taken  at  1  kHz  from  manufacturer's  data  sheets  I 
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Amplifier  Noise  Specifications 


;3.  Devices  shown  in  italics  are  discrete  transistors;  all  other  devices  are  integrated  circuits. 

■4.  uPwr  stands  for  micropower  (very  low  power  consumption)  |  ~J 

'5.  Mbs  devices  have  large  1/f-noise  components  in  the  voltage  noise.  The  1  kHz  values 

! 

cited  are  still  on  the  1/f  portion  of  the  curve.  The  1/f  portion  of  the  voltage  noise 

: 

curves  for  the  other  devices  are  below  (and,  in  some  cases,  well  below)  1  kHz. 

|6.  The  LTC1052  (MOS)  device  is  a  chopper-stabilized  amplifier  (roughly  speaking. 

the  input  is  chopped  at  a  high  frequency,  then  amplified,  then  detected 

synchronously  with  the  chopping  frequency).  It  is  representative  of  the 

best  low-frequency  performance  obtainable  with  a  MOS  device.  Essentially, 

- — - 

there  is  no  1/f  region  in  the  voltage  noise.  j 

|7.  All  of  these  devices  require  biasing  and/or  feedback  resistors.  The  effects  of  these 

I 

resistors  must  be  considered  when  determining  the  overall  amplifier  noise. 

] 

In  most  cases,  low-current-noise  amplifiers  can  be  designed  for  minimal 

I 

effect  from  biasing/feedback  resistors;  for  very-low-voltage-noise  amplifiers. 

the  noise  from  those  resistors  can  limit  the  achievable  noise  performance. 

1 8.  The  two  BJT  discrete  devices  are  shown  for  two  operating  points  each.  One  of  the  1 

I 

advantages  of  using  discrete  devices  is  that  the  noise  performance  can  be 

. i _  . 

tuned  by  means  of  the  collector  (or  drain)  current.  Throughout  this  adjustment 

...j _ 

range,  the  eAiA  product  is  roughly  constant.  | 

1 9.  Very  low  voltage  noise  is  generally  obtained  by  massively  parallel  emitter  reaions  in 

: 

BJT  devices.  This  has  two  important  consequences:  (1)  the  input  capacitance 

- - 

is  large,  and  (2)  the  devices  are  relatively  large.  The  second  point  is  not 

important  unless  the  amplifier  must  be  integrated  onto  a  chip  with  limited 

real  estate:  for  example,  a  single  MAT02  occupies  2x2  mm. 

10,  The  "noise  impedance"‘eA/iA  is  given  for  each  device  but  should  not  be  misinterDreted.  I 

This  number  is  often  cited  as  the  resistance  to  which  the  source  resistance 

must  be  matched  for  best  noise  performance.  It  is  not,  however,  good  practice  • 

=  - 

to  select  an  amplifier  on  the  basis  of  the  eA/iA-to-source-resistance  match. 

If  the  source  (transducer,  for  example)  fs  primarily  resistive,  then  the  amplifier's 

noise  contribution  is  negligible  compared  to  the  Johnson  noise  of  the  source 

over  a  range  of  resistance  centered  on  the  eA/iA  value.  The  width  of  the  region 

to  either  side  of  eA/iA  is  given  (roughly)  by  the  eAiA/4kT  value.  If  eAiA/4kT 

is  0.01  and  eA/iA  is  100  000  ohms,  then  the  amplifier's  contribution  is  less  than 

the  source's  Johnson  noise  for  source  resistances  of  1000  ohms  to  10  megohms. 

It  is  never  appropriate  to  "match"  amplifier  noise  resistance  to  the  magnitude 

of  a  source's  impedance  if  that  impedance  is  primarily  reactive.  I 

1 1 .  The  quantity  eAiA/4kT  is  the  product  of  noise  voltage  and  noise  current  normalized 

by  4  times  Boltzmann's  constant  times  absolute  temperature.  The  product 

4kT  is  approximately  16  x  10''21  in  SI  units  at  room  temperature.  I 

— 

12.  Current  noise  is  exponentially  dependent  on  operating  temperature  in  JFET  devices. 

The  values  given  are  generally  at  25  or  30  deg.  C.  Significant  degradation 

n  current-noise  performance  can  be  expected  with  operation  at  elevated 

temperatures. 

L 
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Sound  Speed  in  an  TH, 


Thermal  Conduction  and  Viscosity  Thermal  (Fourier)  Diffusion  Equation 
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its  normal  in  the  y-direction 
=  Shear  viscosity  [kg/m-sec) 
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Diffusion  Equations  Evanescen 
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Kinetic  energy  density  =  (1/2)  p„u’ 
Positive-definite  energy  dissipation 


Thermoacoustics  is  superficial  science 
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Thermoacoustic  Refrigerators  DEfTAF. 


DeltaE  Model  Agreement  Thermoacoustic  Refrigerators 
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WHAT  IS  NONLINEAR  ACOUSTICS? 
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FIGURE  6-9 

NOISE  PULSE  I  AT  VARIOUS  DISTANCES 
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O'nit)  ~ 

bnit)  =  slowly  varying  function 
=  constant  in  linear  theory 
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Particle  Displacement  Wavefields  NORMAL  MODE  EXPANSION 


HAMILTONIAN  FORMALISM  POTENTIAL  ENERGY 
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Theory  based  on  5  isotropic  constants:  p,  K,  A,  B,  C 
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EVOLUTION  EQUATION 
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x=18.3  mm 


solid  lines:  experiment 
dashed  lines:  theory 
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CHAOS  Contents 
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Physicist:  Determinism 
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Acoustics 


Vihral.c’d  iKniid  lay('r 


after  Ch.  Merkwirili,  Gottingen 
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Vibrated  Ikniid  layer 
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Basic  notions:  state  space  Basic  notions:  attractors 


Example:  resonance  curves  turning  over, 
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Basic  notions:  strange  attractor  Basic  notions:  Bifurcation 


Basic  notions:  Bifurcation  Basic  notions:  Bifurcation  diagram 
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Frequency 


Basic  notions:  Lyapunov  Exponents  Basic  notions:  Lyapunov  Exponents 
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give  underdetermined  results? 
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Relation 

Theory  Experiment 


Nonlinear  time  series  analysis  Nonlinear  time  series  analysis 

Embedding: 


155 


Chaotic  attractors  usually  have 
a  fractal  dimension,  e.  g.  D  =  2.35 
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Nonlinear  Time  Series  Analysis  Nonlinear  Time  Series  Analysis 

Dimension  Estimation 


Nonlinear  time  series  analysis 
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Nonlinear  Bubble  Dynamics  Nonlinear  Bubble  Dynamics 

Bubble  Model  Parameters  Bubble  Model 


Nonlinear  Bubble  Dynamics  Nonlinear  Bubble  Dynamics 

Radius  -  Time  Curves 

Resonance  Curves 

excitation  pressure 
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Period  1  Solution  from  Main  Resonance 
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Bubble  Radius  at  Rest:  55  /iin,  3/1  resonance, 
upper  branch,  period  doubled 

Coexisting  Solutions  from  2/1  and  3/1 
Resonance 


Rn[nm] 

Intermediate  High  Driving 
Small  Bubbles 
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Nonlinear  Bubble  Dynamics  Nonlinear  Bubble  Dynamics 

Radius  -  Time  Curve,  20  kHz,  80  kPa 


Nonlinear  Bubble  Dynamics  Nonlinear  Bubble  Dynamics 

Radius  -  Time  Curves,  20  kHz  Radius  -  Time  Curves,  20  kHz 

Rn  =  8.95  urn  Pa  =  90kPa  (20/1)  Rn  =  7.44  urn  Pa=115kPa  (20/1) 


uu/u 
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Bubble  Radius  at  Rest:  8.77  tim,  100  kPa  Bubble  Radius  at  Rest:  6.71  pm,  130  kPa 

Solutions  from  the  20/1  Resonance  Solutions  from  the  20/1  Resonance 


High  Driving,  Small  Bubbles 


Nonlinear  Bubble  Dynamics 
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Nonlinear  Bubble  Dynamics 


Single  Bubble  Sonoluiniii(‘S(  ('uc(' 
Spherical  flask 
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Single  Bubble  Sonoluminesceiice  Bubble  Sonoluminescence 

Square  container  Cylindrical  container 
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Laser-produced  bubble 

in  water  Acoustic  cavitation 


I 

pci 


171 


Frequency 


172 


[TR-66] 


174 


1 


Quantum  Mechanics  Tutorial 

J.  D.  Maynard 

The  Pennsylvania  State  University 

Purpose 

How  are  Quantum  Mechanics  and  Acoustics  related? 

•  Both  involve  the  same  wave  phenomena  (Wave  Mechanics) 

•  The  attenuation  of  sound  involves  molecular  collisions  which  must  be 
treated  with  quantum  mechanics 

•  Solid  crystals  have  quantized  sound  waves  (Phonons) 

•  Interaction  between  phonons  and  electrons,  effects  of  pressure,  magnetic 
fields,  light  (Raman,  Brillouin  scattering),  etc. 

•  Acoustics  in  Macroscopic  Quantum  Systems  (Superfluids) 


2 


History 


Events  which  led 


to  Quantum 


Mechanics 


•  Plank's  theory  of  Blackbody  Radiation  (1901) 

-  Due  to  quantized  modes  of  a  solid  (phonons) 

-  Not  due  light  behaving  as  particles  (many  textbooks  imply  this;  Plank 
did  it  right) 

-  Energy  of  mode  is  quantized  as  Plank’s  constant  (h  =  6.625  x  lO"^'* 
J-s)  times  frequency  of  mode 


176 
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Events  which  led  to  Quantum  Mechanics,  continued 

•  The  Photoelectric  Effect  (1905):  Light  causes  electrons  to  be  emitted 
from  a  metal 

—  Emssion  begins  almost  immediately,  even  for  light  intensity  of  only 
10-iow/m2 

—  The  energy  of  the  electron  is  proportional  to  the  light  frequency  v 

—  Einstein:  Light  waves  act  like  particles  (photons)  with  energy  hv 

—  Many  textbooks  say  photoelectric  effect  shows  that  light  must  act 
like  particles:  absolutely  incorrect  [Scully,  Phys.  Today,  Mar  1972] 

-  Although  wrong,  lightwaves  acting  like  particles  historically  suggested 
that  particles  might  act  like  waves 


intraviolet  light 
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Events  which  led  to  Quantum  Mechanics,  continued 

•  Rutherford  Scattering  (1911):  Electrons  orbit  heavy  nucleus 

•  Bohr’s  Theory  of  Atomic  Spectra  (1913):  atoms  emit  light  at  discrete 
frequencies 

—  Electrons  orbit  nucleus  classically,  except  only  at  certain  radii 

-  The  allowed  radii  are  such  that  the  electron’s  angular  momentum  is 
an  integer  multiple  of  Plank's  constant,  divided  by  27r. 
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Events  which  led  to  Quantum  Mechanics,  continued 

•  Compton  Scattering  (1923):  wavelength  of  light  is  shifted  when  scat¬ 
tered  from  an  electron 

—  Historically,  gave  further  evidence  for  particle  nature  of  light 
—  Like  photoelectric  effect,  conclusion  is  absolutely  incorrect 

•  De  Broglie  Waves  (Ph.D.  Thesis,  1924):  Electron  (momentum  p,  en¬ 
ergy  E)  is  a  wave  with  wavelength  A  =  h/p  and  frequency  v  =  B/h. 
Explains  all  preceding  experiments,  and  predicts  that  electrons  diffract 
like  waves 

•  Schrodinger’s  Wave  Equation  for  particles  (1925)  [more  later] 

•  Heisenberg's  Uncertainty  Principle  (1927)  [more  later] 

•  Davisson  and  Germer  (1927):  Electrons  in  a  crystal  lattice  diffract 


Myths  and  the  Mystique  of  Quantum  Mechanics 


•  Wave-Particle  Duality  :  Electrons  may  behave  either  as  waves  or  par¬ 
ticles 

—  The  Law  of  Physics  is  the  Schrodinger  Wave  Equation;  everything 
must  be  explained  in  terms  of  waves 

—  The  Classical  Particle  picture  is  only  an  approximation  (e.g.  ray  trac¬ 
ing  in  optics) 

—  For  some  reason,  people  are  reluctant  to  give  up  the  notion  of  parti¬ 
cles 


•  Paradoxes  arise  from  Quantum  Mechanics:  There  are  no  paradoxes; 
they  only  arise  if 

-  One  insists  on  giving  objects  particle-like  attributes 

One  uses  “detectors"  which  do  not  obey  the  laws  of  physics 

•  The  Uncertainty  Principle  :  This  has  significance  only  if  one  refuses  to 
give  up  the  notions  of  particles 

•  Quantum  Chaos  :  The  Schrodinger  equation  is  linear;  there  is  no  chaos 
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A  Valid  Picture  of  Quantum  Mechanics 

1.  Solve  a  Boundary  Value  Problem  which  is  mathematical,  rigorous,  and 

has  a  unique  solution 

I'JOTE:  After  step  1)  absolutely  nothing  happens! 

2.  Make  a  Measurement 

•  A  measurement  involves  a  very  large  number  of  complicated  elements 

•  The  minimum,  simplest  element  is  a  graduate  student  with  lab  note¬ 
book 

•  A  system  with  many  elements  is  sensitive  to  small  perturbations 

•  No  system  is  totally  isolated  (no  shield  for  gravity  waves);  there  are 
always  small  perturbations 

•  The  best  one  can  do  is  use  the  results  of  1)  to  calculate  the  probable 
outcome  of  a  measurement  (  John  von  Neumann  ) 

•  Reference:  Zurek,  Phys.  Today,  Oct.  1991,  ‘‘Decoherence  and  the 
Transition  from  Quantum  to  Classical" 
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Setup  for  Quantum  Mechanics 


Classical  Mechanics 


Newton’s  Law 
Linear  Momentum 
Angular  Momentum 

Kinetic  Energy 

Potential  Energy 
Conservation  of  Energy 


F  =  rtid  —  dpidt 


p  =  mv 
L  =  fx  p 

T  =  —mv^ 
2 


V 


=  / 


F-df 


T  +  V  —  Constant 
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Classical  Mechanics,  continued 


Electrodynamics 

Scalar  and  Vector  Potential 


E<&'M  Force  Potential 


Classical  Mechanics,  continued 

Generalized  Coordinates  n  =  fi  (91  •  •  •  gN,t) 

fw  =  rM  (gi  -  •  ■  97V) 

Because  of  constraints,  N  <  3M 
Velocity  Dependent  Potentials  (91  •  •  •  g/v,  91  •  •  •  97v) 

Lagrangian  L  =  T  -  U  =  L(gi,gi,t) 

Lagrange's  Equations  —  ( — ^  -  —  —  0 

dt  \dgij  dgi 

Note  that  the  first  term  is  second  order. 
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Classical  Mechanics,  continued 


Generalized  Momenta 
Example:  Charge  in  magnetic  field 

The  Hamiltonian  H {qi,pi,t)  =  -  L 

3 

Example:  Charge  in  magnetic  field  H  =  (p- 

2m  V  c  / 


dL 


q 

Px  =  mx  H — Ax 
c 


Hamilton’s  Equations 


dff 


NOTE:  There  are  twice  as  many  equations,  but  they  are  first  order. 
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Classical  Mechanics,  continued 


Conservation  Laws:  Suppose,  perhaps  because  of  some  symmetry,  the 
Hamiltonian  does  not  depend  on  some  qj.  Then 

dff  ^ 

Pj  =  =  0  pj  z=  constant 

In  particular,  if  L  dosen’t  depend  explicitly  on  time,  then  ff  =  constant. 

the  transformtions  defining  the  generalized  coordinates  do  not  depend  ex¬ 
plicitly  on  time,  and  the  potential  energy  is  velocity  independent  (conservative 
forces)  then 


ff  =  T-\rV  =  Total  Energy 
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Vector  Spaces  and  Function  Spaces 


A  =  HiAiXi,  Inner  Product;  A-B  =  HiAiBi 


Functions 

Inner  Product 

Orthogonality 
Unit  Vectors 

Operators 

Matrix  Representation 


i’A  (x)  =  e  5®*  (81^  —  12x) 


<^a\'>Pb  >=  j  “Ip A  (x)  xpB  (x)  p  (x)  dx 

<  V’A  I  t/’b  >=  0 
Xpi  I  %f)j  Sij 

Q 

\'^A>=  P\‘^A>  Example:  -P  =  ^ 
Pij  =<-$i\P\  ■ipj  > 


Expectation  Value 


<  P  >  =  <  i}A\P\^A> 


Commutator 


[P,Q]  =  PQ-QP 


Quantum  Mechanics 


Formal  Theory 


•  At  a  time  to,  the  state  (given  a  label  "^4")  of  a  quantum  mechanical 
system  is  given  by  a  point  in  a  linear  vector  space:  ipA(to)- 

•  The  state  is  normalized  at  all  times:  <  ipA  (*)  I  i’A  (i)  >=  1- 

•  There  exists  a  time  translation  generator  (linear  operator)  T(t,to)  such 
that 


‘tpA  (t)  =  T  (t,  to)  “tpA  (to) 

•  Define  an  operator  H  (t)  such  that  T  (t  +  dt,  t)  =  1  —  iH  (t)  dt/H. 
Since  ^  (t  +  dt)  =  ^  (t)  -}-  {dip/dt)  dt,  then 


H(t)rPit)  =  iTi-^ 
at 
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Formal  Theory  of  Quantum  Mechanics,  continued 
Consider  the  dynamics  of  the  expectation  value  of  an  operator  P : 

=<<Ht)  I  >  -  <  £^(()  I  p^(()  >=  4  <  IP.H]  > 

^  tfi  in  ih 

•  The  Correspondence  Principle 

>=  where  Pc  is  the  classical  quantity 


If  P  =  g, 


ajfc 

dpi 


If  P=.p, 


dHc 

dqi 


These  conditions  can  be  satisfied  if  we  let  H  be  the  classical  Hamiltonian, 
but  with  the  qi  and  pi  replaced  with  operators  satisfying 


fe,gi]  =  0,  b«,Pi]  =  0,  but  [qi,Pj]=iTiSij 


Formal  Theory  of  Quantum  Mechanics,  continued 
There  is  more  than  one  way  to  define  the  g,  and  p,  operators. 


Coordinate  Representation 

9i  qi, 

-r  9 

Pi  -^n— 

dqi 

Momentum  Representation 

Pi  ->  Pi, 

•X  ^ 

Qi 

Opi 

Heisenberg  Representation 

Qi,  Pi  -> 

Matrices 

Example:  [x,p,]  f  (x)  =  x  /  (x)  -  xf  (x) 

=  —xiTi-^  +  ihf  (x)  4*  ihx-^f-  =  iHf  (x) 
dx  dx  ^  ' 

In  3-D  coordinate  representation:  p^l2m  =  — 

2m 


183 


17 


Formal  Theory  of  Quantum  Mechanics,  continued 

©2 

Suppose  H  =  T-\-V  =  —-\-V{r,t) 

2m 

Then  we  get  the  Schrodinger  Equation: 

Hvl/  (f,  t)  =  (f,  t)  +  V  (f,  t)  M/  (r,  t)  =  a— 

2m  at 

Suppose  V  (f,t)  =  V(f).  Separate  variables:  M/  (f,t)  e~^K 


Then  we  get  the  Time-independent  Schrodinger  Equation: 


-  —  (fo y  (0  V*  (r)  =  ■EV' (0 

2m 


where  E  =  ^cv.  For  closed  systems  E  will  be  a  discrete  eigenvalue,  or  Quan¬ 
tized  Energy  Level.  The  eigenfunction  ^(f)  is  called  a  Stationary  State. 
The  eigenfunction  for  the  lowest  E  is  called  the  Ground  State.  Discrete 
eigenvalues  are  labeled  with  integers  called  Quantum  Numbers. 
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Examples 


Free  particle:  1/^  =  0 


One-Dimensional  Free  Particle 


2m  dx^ 


=  Eip,  E  =  Tiuj 


-fA:V  =  0, 


dx^ 


k  =  —y/2mE 
Ti 


Note  classical  momentum  p  =  V2mE  -)•  p  =  Tik. 


Solutions  (x,  t)  = 

Waves  of  wavelength  A  =  2Tr/k  =  h/p  and  frequency  u  =  a;/27r  =  E/h  (de 
Broglie!). 


Note:  cv  =  E/H  =  {'h/2m)  k^  =  u  (k)  Dispersion 

For  a  common  acoustics  wave,  w  =  cok.  Linear  Dispersion 

Dispersion  is  a  fundamental  difference  between  free  particle  Schrodinger  waves 
and  common  acoustic  waves. 
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Wave  Velocities  and  Dispersion 

For  v|/  (x,t)  =  (kx  -  ojt)  =  </>  is  called  the  phase. 

What  is  the  relationship  between  x  and  t  such  the  phase  ^  is  constant? 

,  dd)  dx 

<f>  =  constant  -4  —  =  0  =  k - u 

dt  dt 

Solving:  {dx/dt)  =w/k  —  Vp^the  Phase  Velocity. 

Suppose  xj/  (x,  t)  -  cos  (kx  -  ut)  +  cos  [(k  +  A/t)  x  -  (w  +  Aw)  t] 

=  2  cos  i  (Akx  -  Awf)  COS  {kx  —  (Dt)  =  envelope  x  carrierwave 


The  envelope  moves  with  velocity  Aw/A/k  duj/dk  =  Vg,  the  Group  Velocity. 
For  linear  dispersion,  Vg  —  Vp.  Otherwise,  vg  ^  Vp. 


Wave  Superposition  and  Wave  Packets 

The  General  Solution  is  a  linear  combination  (superposition)  of  eigenfunc¬ 
tions 


A  (k) 


Suppose  A(k)  is  peaked  around 
ko,  as  shown  in  the  figure: 


Note  that  xj/(x,<  =  o)  is  the  Fourier  Transform  of  A(k),  which  is  a  wave 
packet: 
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Wave  Packets,  continued 

From  a  theorem  for  Fourier  Transforms,  the  product  of  the  width  of  the  curve 
in  k  -space  and  the  width  of  the  envelope  in  x  -space  is  a  constant.  If  the 
curve  in  k  -space  is  a  Gaussian,  exp  [-  (fc  -  fco)^/2] ,  then  the  envelope  in  x 
-space  is  also  a  Gaussian,  and  the  product  of  the  widths  is  minimized. 

Recall  that  A(k)  is  peaked  at  ko  Let  K  =  k-kQ  Then 

VI/  (x,  <  =  0)  =  j  A  (kye^'^^dk  =  J  A  (ko  +  k)  e^'^dn 

For  times  f  >  0,  we  approximate  u  (k)  near  A:o  with  a;(fc)  ~  wo  +  (du>/dk)  k  = 
4-  VgK.  Then 

vi^  (x,  t)  =  J  A  (ko  +  k) 

^  ei(k,v,-^,)t^ik^(z-v,t)  J  A(ko  +  k) 

^  [(a;  _  ^  0] 

Taking  the  modulus  eliminates  the  phase  factor,  leaving  the  original  envelope 
of  the  wave  packet  evaluated  at  (x  —  Vgt) ,  i.e.  moving  with  the  group  velocity. 
With  non-linear  dispersion,  the  wavenumbers  higher  and  lower  than  ko  move 
with  different  velocities,  causing  the  wave  packet  to  spread. 
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Examples 


Piecewise  Constant  Potentials 


Piecewise  constant  potential:  V  (x)  =  Vi  =  constant  for  i,  <  i  <  Xj+i 


For  Xi  <  X  <  Xj+i, 


+  W  =  E  =  7iu; 

^  +  kfip  =  0,  ki  =  \s/^m{E-Vi) 


Solutions:  Ur  (x,  t)  = 


Note  that  in  some  regions,  E  may  be  less  than  V5,  and  ki  will  be  imaginary. 
Solutions  of  the  form  -^{x)  oce“'“  are  called  Evanescent  Waves. 


The  coefficients  Ai  and  Bi  are  found  by  satisfying  the  conditions  that  V’(®) 
and  its  derivative  d-tl^/dx  be  continuous  at  boundaries. 

Satisfying  all  boundary  conditions  may  result  in  solutions  existing  only  for 
discrete  values  E,  which  are  the  quantized  energy  levels. 
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Examples 


step  Potential 


E- 

V(x)  =  0- 


-V(x)  =  Vo 


X  =S  0 


On  the  left:  ^  (x)  =  e‘^  +  Re  k  =  —V2mE 

Ti 

On  the  right:  V  (®)  =  Te’’®,  q  =  4\/2m  {E  —  Vb) 

Ti 


Equating  the  V’  and  the  d-^/dx  from  the  left  and  right  gives  two  equations, 
which  are  solved  for  the  two  unknowns  R  and  T,  the  Complex  Reflection 
and  Transmission  Coefficients. 

R={k-q)/{k  +  q)  f^2k/(k  +  q) 


If  E  >Vo,  then  there  are  waves  on  both  sides,  but  with  different  wavelengths. 
If  a  wavepacket  were  incident,  part  would  be  reflected,  and  part  transmitted. 
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step  Potential,  continued 


E 

V(x)  =  0 


x  =  0 


V(x)^Vq 

E 


If  E  <Vo,  then  g  -4  iq,  and  the  solution  on  the  right  becomes  an  exponentially 
decreasing  evanescent  wave. 


The  R  and  jT  become  more  complicated  complex  numbers,  but  |  R  p,  the 
energy  reflection  coefficient,  would  be  unity. 

If  a  wavepacket  were  incident,  the  superimposed  waves  with  E  >Vo  would  be 
completely  reflected. 


X  K  0  X  »  o 


On  the  left:  ^'(x)  =  e**“  +  k  =  V2mB/fi 

In  the  middle:  rp  (x)  =  g  =  y/Tm^E^^^/Ti 

On  the  right:  V’ 


[f  o  = 


9  =  *«, 


— i  sinh  (ga) 


2/:/ccosh  (ga)  -i(k'^  -  k?)  sinh  (gaV 


T=  etc. 
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Examples 


Square  Well  Potential,  E  <  Vo  (Bound  States) 


$ 

€) 

& 


0-i- 

- 1 - ^ - 1 - 

^3C*s  —  eft  **0  JTMIxft 

On  the  left; 

^  (x)  =  Ae’^,  K  =  y/2m  (Vo  -  E)/7i 

In  the  middle: 

=  Scoskx-i-Cs\nkx,  k  =  V2TnB/Ti 

On  the  right: 

^  (x)  =  De~'^ 

The  boundary  conditions  require  either  C  =  0  and  ifctan  {ka/2)  =  k,  or  B  =  0 
and  kcot{ka/2)  =  —k.  These  two  possibilities  correspond  to  symmetric  and 
anti-symmetric  bound  states. 

The  transcendental  equations  yield  a  finite  number  of  discrete  values  for  the 
energy  levels  E. 

If  Vb  oo,  then  the  eigenfunction  V"  must  vanish  at  x  =  ±a/2  (the  derivative 
dxl)/dx  will  be  discontinuous).  The  eigenfunctions  are  the  same  as  a  string 
clamped  at  the  ends,  with  fc  oc  n,  an  Integer.  However,  the  string  natural 
frequencies  are  proportional  to  n,  whereas  the  quantum  energy  levels  are 
proportional  to  n^. 
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Square  Well  Potential,  E  <  Vo  (Bound  States),  continued 


n 

s. 

u 

A 

f\  1 

y 

H 

■a 

For  finite  Vb,  the  eigenfunctions  are  shown  with  dashed  lines,  and  the  energy 
levels  are  shown  in  (a). 

For  infinite  Vb,  the  eigenfunctions  are  shown  with  solid  lines,  and  the  energy 
levels  are  shown  in  (b). 


193 


37 


Examples 


The  Simple  Harmonic  Oscillator  Potential 


The  Simple  Harmonic  Oscillator  Potential:  V  (x)  = 

where  K  is  the  spring  constant,  m  is  the  mass,  and  u»o  is  the  classical  oscilla¬ 
tor’s  natural  frequency. 


The  Schrodinger  Eq. 


Let  y  =  y/mi^/n  I  and  let  V  (i)  =  exp  (-y^/2)  k  (y) .  Then 


d^h  dh  ^  /  E 

dy^ 


h{y)  =  0 


If  a  series  solution  for  h  {y)  is  tried,  then  the  series  diverges  unless  (E/Tujjo  -  1/2)  = 
n,  an  integer.  Thus  the  energy  levels  are  quantized  with 

En  =  ^  Tuoq 

The  eigenfunctions  hn  (y)  are  the  Hermite  Polynomials. 
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The  Quantum  Mechanical  Simple  Harmonic  Oscillator,  continued 
Normalized  Wavefunctions: 

vi/„  (x,  t)  =  x)  e-(W7i)-’/2 


'BHBI 

(y)  =  2y  A _ 

_ 3Hlj/2 

InHHIHl 
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The  Quantum  Mechanical  Simple  Harmonic  Oscillator,  continued 

For  a  classical  simple  harmonic  oscillator,  x(t)  =  ^  cos  (wot),  and  the  total 
energy  \s  E  =  (1/2) 

Classical  result:  y  (t)  =  yJlE/nujl  cos  (wot) 

Quantum  Expectation  Value  for  state  M/„(i,t)  :  <  1 1 1  >=  o 


Measurement  couples  states:  M/  (i.f)  =  \|/„  (x,t)  +  (x.t) 

Now  <  M/  1 1 1  \|/  >=<  1 1 1  x|/„+i  >  +  <  y/„+i  I  X  I  M/„  > 

=  2Re  <Wn\x\  > 

=  {i’n  (x)  (v-n+i  (x)  e-’("+t)‘^‘)  dx 

=  2Re  j  ipn  (x)  xtfjn+i  (x)  dxe“^‘ 

=  2Rey/ nTif2rriu}Qe^^^^^ 

=  J2nTuuo/mcjQ  cos  (wot)  ~  -J 2E{m(jjQ  COS  (wot) 


Conservative  Central  Force:  F  =  Ff  ->  V  (r)  =  V  (r) 


Classical  consequences  for  angular  momentum,  L  =  rxp. 


dL  ^  ^  ^ 

——  =  T  =  rxT=0 
dt 


L  =  constant 
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Central  Force,  continued 


Quantum  Mechanics: 


Hi})  =  — - — V”  (r)  V"  =  Eij; 
2m 


_  7l2  1  5  /  \  1  d  ^d\  1  02  1 

2mr^dr  dr)  ^  2mr^  [sin 050  V  " ^50 )  sin^^J  +  ^  ^ 


Separate  variables:  ^{r,e,<f>)  =  (l/r)  ii  (r)  ©  (0)  <t>  (^) .  Plug-in  and  divide: 
1  (P<x> 

—  constant  =  —mr  O  (^)  = 

Single  valued  V-  (r,  0,  -f  2ir)  =  V-  (r,  0,  m  =  0,  ±1,  ±2,  ±3  -  - 

1  d  f  ,  ^dQ\ 

©^^  i"'"  =  constant  =  -/(/-!-  1) 

Solution:  Legendre  Polynomial,  ©(0)  =  f^(cos0)  and  Z  =  integer  <|  m  | . 


d^jR  rz(Z  -f- 1)  I 

“5::^+  -SJ 


Central  Force,  continued 

Combine  0,  <^  ;  y;-  (0,  (-1)™  (cos  0) 

Quantum  Angular  Momentum  L  =  r  x  p  =  -i^r  x  V  =  Lix  -}-  LyyLzz. 


L. - a  (y|  - 


•]  =  (  sin^-^-f-cot0cos<^— 
/  \  Op  dd 


L,  =  -<7i  -  x|)  =  iS  (-cos.^1  +  cot8 sin 


„  Tl^  1  d  f  ^dih\  l2 


Note  —  0,  [H,Lt]  =  0  >  z=  const,  and  <  L,  >  =  const. 

Note  L2j.m  (0^  ^  ^  30^  ^  (0,  <^) 

Expected  modulus  of  Angular  Momentum:  L  =  V<  >  =  •y/r(Z'+T)7i 
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Example 


Coulomb  Central  Force:  V  (r)  =  (Ze)  ejr. 


2m 


2mr2 


+  V(r) 


R=-ER 


—R(r)  =  e  {2knr)^Gn  (2fc„r) ,  Laguerre  Polynomials 

T 

1  _  mZ^e^ 

kn  =  -y/7^n,  B,,  =  -^,n=  1,2,3,- ■■ 

n  2Ji 

For  a  given  n  :  Z  =  0, 1,2  •  •  •  (n  -  1) ;  m  =  -Z,-Z  +  1,- • -Z  -  1,Z,  [  (2Z  +  1) 
values] 

Historical  nomenclature:  Z  =  0. 1,2,3  are  referred  to  as  ‘‘s,  p,  d,  f"  states. 


Energy  eigenvalue  £?„  depends  only  on  n.  Eigenfunctions  for  different  Z,m  are 
^degenerate.  Number  of  degenerate  states  =  (21  +  1)  =  n^. 


Arbitrary  z-direction  not  a  problem,  because  sum  over  degenerate  states  is 
isotropic. 
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Quantum  Mechanics 


History  and  Formal  Theory,  continued 


•  Stern-Gerlach  Experiment  (1922)  Electron  in  a  inhomogeneous  mag¬ 
netic  field 


Deflection  magnetic  moment  -4-  charge  with  angular  momentum.  But 
atoms  in  ground  state  had  I  =  0  electron's  orbital  angular  momentum  = 
0! 


•  Dirac  (1927)  Extended  formal  Hamiltonian  theory  to  include  relativistic 
dynamics  -4  Electron  has  intrinsic  angular  momentum,  Spin 

Spin  angular  momentum:  |  5  |=  y/s  (s  -f  l)7i,  with  s  =  1/2.  Note  2s-f  1  =  2. 
Z-component  of  electron  spin:  Sz  =  m,7i,  m,  =  -^  or  -f  i 

Now:  Quantum  numbers  for  the  single  electron  atom:  n,l,m,ms 

Wave  function  gets  a  two-element  vector  (  Spinor  )  attached.  Spin  operators 
are  2x2  matrices  (linear  combination  of  four  Pauli  spin  matrices  ). 
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Quantum  Mechanics 


Formal  Theory,  continued 


For  systems  with  many  particles,  there  are  generalized  coordinates  (and  op¬ 
erators)  for  each  particle. 


For  classical  identical  particles,  the  particles  may  be  distinguished  by  following 
their  (precise)  classical  trajectories. 


In  quantum  mechanics,  particles  are  waves,  and  the  Principle  of  Superposition 
means  that  identical  quantum  particles  are  Indistinguishable 


Note:  Any  complete  set  of  dynamical  variables  (or  quantum  numbers)  K 
which  describes  a  single  particle  can  also  be  employed  for  n  particles  of  the 
same  kind,  even  if  the  particles  are  interacting. 


Consequence:  vi/  (//particles)  involves  products  |  ifi  >|  K'2  >  •  •  •  |  Kn  > 

Since  the  particles  are  indistinguishable,  exchanging  any  pair  of  |  fC  > 's  must 
give  the  same  result  (  Exchange  Degeneracy  ).  There  are  only  two  forms 
for  XU  for  which  ]  vu  p  is  invariate  under  any  exchange.  Let  P  indicate  one  of 
Nl  permutations.  Then: 


VUs  =  I  >1  K2  >  •  •  •  1  > 

P 


^  ap  I  /Cl  >1  i<'2  >  •  •  •  I  Kn  > 

p 

where  ap  =  1  for  even  permutations,  and  =  -1  for  odd  permutations. 
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Formal  Theory,  Identical  Particles,  continued 
Note:  Hfs  (-Ki-  Kj .  •)  =  (•  •  Kj  ■  Ki  •  •) 

and  Ml  A  (•  •  Ki  ■  -Kj  ■  ■)  =  -vl/^  (■  •  Kj  ■  -Ki  ■  •)  . 

Note:  =  0  if  any  two  particles  have  the  same  Ki. 

•  History:  Pauli  Exclusion  Principle  (1925).  [From  experimental  obser¬ 
vations  on  multi-electron  atoms]  There  can  never  be  one  electron  in  the 
same  quantum  state. 

Consequence:  Electrons  must  have  anti-symmetric  wave  functions, 

Generalization:  Particles  with  half-integral  spin  quantum  numbers  must  have 
antisymmetric  wave  functions;  such  particles  are  called  Fermions  .  Particles 
with  integral  spin  quantum  numbers  must  have  symmetric  wave  functions; 
such  particles  are  called  Bosons  . 

Thermodynamic  Distribution  functions: 

Fermi  Statistics  (B)  =  ^ 

Bose  Statistics  *  (E)  =  ^ 
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Perturbation  Theory]  Time  Independent  Potential  Fields,  V (r) 


Suppose  one  can  solve  /foV’n  =  •®n'0n 

Wish  to  solve  (Ho  V)  ^„  =  E„V»n-  For  small  V  : 

cr  e;+ <  I V I  <  >  +  2  + . . . 


If  degeneracy  then  diagonalize  matrix. 


Variational  Method: 

=  Eoipo  is  equivalent  to  finding  the  ^  which  minimizes 

<‘ip  \ip  > 

Write  a  -ip  with  parameters,  and  minimize  SH  with  respect  to  the  parameters. 
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Perturbation  Xheory  Time  Evolution  (V  may  depend  on  time) 


Wish  to  soive  {Ho  +  V)  M/  =  i7i  {d’^/dt') 

If  W°  =  7ia;„V'2.  then  \J/0  (f)  = 

n 

Assume  (*)  =  ^  Ai  (*) 

n 

Piug  in:  ^  ^  V  <  V-?  K  I  V’S  > 

at  tfi 

n 

Assume  c,  (f  =  -oo)  =  1,  and  c*  (f  =  -oo)  =  0  : 

(0  =  4  /■  <  *“  I  V  I  V'S  > 

Example:  EM  Radiation  A=  /  A(a;)e*(^'' 

^  — OO 

H  =  i  (p-_  =  ^  -  Xa  .  p  +  |5:v .  A-+ 

2m  V  c  /  2m  me  2mc  2mc^ 

2m  me 

a;  =  wji:  —  u»5  — ^  Resonance  .  <  ^°  |  1  >=  0  Selection  Rules 
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[TR-2I 


Tuniif>^-up  g  Quasi  cry  5  to,  I 


23 

Problem;  Sc/ve  ciccthjn  TffiiH  lO  /aMS 


®  /\nderSoKi  f»caliz«it ion  ®  Universal  cdn4i4ctAnce  fluctiutioMS 

o  Khrenov '  BoKm  effect  ®  Mormnl  clettrM  persis't’cnt  carren'ts  • 


[TR-3] 


“  Phase  coherence  dii  the  rc«le  of  microns 

Megascopic  Phase  cdtereace  on  the  scale  of  miillanj  af  microns 

Eripenmest*  :  •  fhme  coh trace  m  a  l-D  wire  10  m  long 

•  Density  of  states  m  a.  i^uasi  crystal  >■  Iw  m  4iaai«ter 

Classical  (acowtic*)  analog  Systems  ("analog  Computers"' ^ 

Advantages  •  Prease  analsgs — V(r)3'P  ~  ® 

®  All  coiwlitioas  and  parameten  may  be  precisely  controlled 
or  measured 

•  Precise  measareoient  of  cigenraluec^  eigenfunctions^ 

density  of  states  ^  etc. 

•  May  study  ftmc' dependent  or  non- linear  effects  exactly 
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[TRs  4  &  5  Unavailable  At  Time  Of  Printing] 

[TR-6] 

R^PMriON  WITH  AHP£RS0N  localuation 
UffOTHETICAt,  rnOPLEMIN  NOISC  RCPKCTION 


[TR-7] 


reKlopic  ARMf  or  IOEHTICAl  mbs  ( E«ise  ol  w<Muf«ct«re) 


£X((^iHpld:  ^(ccfpdiC  a.  rtet^lllc  cr^yj-fal 

7—*'  ® _ ^  ©  ©  ©  @  ©  © 

^  "t  Potifwd  ion  lattice 

EtectricJ  Condactiffily  oC  dfslan^c  bcTwcc/i  scattcrers 


206 


[TR-8] 

Of«.tVE  FROM  frXOUP  THfrtRV; 

Mfllhemutics  -  FleaUETS  THEOREM 
5oliJ  state  -  BtOcH's  THE0^EM 

For  4  iysittn  with  a  p<rioAic 
potent  kI  or  iriptJance  ,  Ihe 
et<|eiifuivct(oii$  art  e>teii<(eJ  : 

where  Urt  + 


IVk  I  jcoflstant  for  aII  X 


Solid  State:  Block  Wave  ^auctions  — ♦■ 


**rl*i| 


[TR-9] 


rewopic  roTCNTiAL  =  v.(xi 

:¥¥3jmHt¥^ 

M 


+  AV  RAN06M 


(b) 

na:utofn  pn(cttll«l  ciicrKjf  intHKluceil  lijr  Andenon;  («)V*  C^)  ^  +  AV 


ANPERSOM  lOC/qH£/iTfQN 


Wave  function  ifi  of  nn  electron  whcii  («)  ICxtcmlcc! 

(A)  localized  3t«te3« 


(  M58) 


rt.  LuI>am  -»  J.  Luscoufte  ^  Fk^f.  f{ey.  E>^5  .  ‘loHS 
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[TR-11] 
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0  MOSAIC  •olufPf  15  Nymbfr  4  W,nt*r  1987/8 


[TR-25] 


EHccts  of 


OlU«t  J’CJ'IOollClty 


CNWCY  CtiO 


spectrum  Density  of  states 


[TR-26  Unavailable  At  Time  Of  Printing] 
[TR-27] 


[TR-29] 


Meo.surew<nt  of  Quat  icry&tal  £iq6nfuntt 


Mil 


Scanner 


[TR-30] 


2 


fTakes  Spotlig^ht 
in  Study  of  Solids 

onu.  - K 


[TR-32] 


PliartKcr  cuJ 


<*iflneArit 


[TR-33] 


[TR-34] 


Does  Nonlinearity  Weaken  iderson  Localization? 


Reference:  ^  ^ 

1.  P.  Devillard  and  B.  Souillard,  J.  Stat.  Phys.  43.  423  (1986) 

Fixed  output  t,  find  t/r  decays  as  power  law  for  strong  nonlinearity  ^ 

2.  B.  Doucot  and  R.  Rammal,  Europhysics  Lett.  3,  969  (1987) 

Fixed  output:  power  law  decay  -  Fixed  input:  exponential  decay  ^ 

3.  C.  Albanese  and  J.  Frohlich,  Commun.  Math.  Phys.  116,  475  (1988) 

Rigorous  theorem:  Eigenstates  of  NLS  eq.  remain  localized  X 


4.  Q.  Li.  C.  M.  Soukoulis  St  Pnevmatikos,  and  E.  N.  Economou.  Phys. 

Rev.  B  11888  (1988)  y 

A  soliton  can  force  its  way  through  a  binary  alloy 

5.  A.  Soffer  and  M.  1.  Weinstein.  Commun.  Math.  Phys. 

Same  as  3.  X 

6.  R.  Bourbonnais  and  R.  Maynard.  Phys.  Rev.  Lett  64.  1397  (1990) 

Superpositions  of  localized  states  spread  due  to  nonlinearity  ^ 

7.  Yu.  S.  Kivshar.  S.  A.  Gredeskul.  A.  Sanchez,  and  L.  Vazauez  Phvs 

Rev.  Lett.  64.  1693  (1990)  ^  ^  y  y 

Same  as  4,  but  only  for  sufficiently  strong  soliton 

8.  R.  Scharf  and  A.  R.  Bishop,  "Nonlinearity  with  Disorder",  cd.  F. 

Abdullaev,  A.  R.  Bishop,  and  S.  Pneuvmatikos  (Springer.  Berlin,  1992)  X  X 

The  nonlinear  Schrodinger  equation  on  a  disordered  chain 
Numerical  re.sulLs;  same  as  7 


217 


t  anAlysii  versus 


(.  <^0(4^1  !(«/  ^ 


\/V  T(co) 

t 

pouirtir  TrwKjj^) 

i  H-Ct) 


»«  ••<  «  « 


jjonlinti^rity  !»  g  SlreichtJ 


^  =  Ma$)  j  liiijfh 


Newton ; 


^rc.  Itd^tk 


*uati\ 


M 

nr^  "  r  =0  w^v? 

T  -  T.  [u  (i)[  .  ,]j 


P f'ffAu.ci  /t^  6  \  / 

Malftces  (»  •  ){*!!^  ,*)  *•**• 


c." 


Us  C;c-tr«  ^<c^rce  4  -Frcc/lfiH^ 

Mtiify  CoViM\OAs  locAiy  ^  within  (s&onJI)  L 


c  .  L,a» 

'  He  n 


WovK  jolJoH  i.  » L. 

He  ff 


Ois+<«it«  tV'a.veleA 


jjgnlme^r  puUo.  .Wt.., 


■j  S^e<A  0+  w»vc  ^  ^  i= 


Fimte  <x*,pl,Wc:  Wa«U, 

U^««uatw„,  SujPcrfU  Hdiuw  UU)  =  \/^  VV^s  ^ 
TliiK'iA  $0«hjI 


[TR-41] 


No  scAtfci^ffrs 


(k 


) 


[TR-42] 


Pd  Ka.titrtr$  «4  fifinjtn^Mr  pul 


s  e 


221 


223 


[TR-49] 
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HTTlE-  TMBdHO^  SLl&t^riy  CM^fLCAJ^ 
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RESONANT  ULTRASOUND  SPECTROSCOPY 


RESONANT  ULTRASOUND  SPECTROSCOPY 
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no.  3.  U)  Summary  of  Um  lowlying  phooon  brtachM  in  U]>.Sr^04  (umplci  MIT-I  ud  MIT* 2).  Tbe  Usm  an  (uidn  U) 
Um  eye.  Tb«  modce  heve  boea  Ubeiod  eooordlni  to  Ref.  4.  The  dUpenioo  curvee  are  oaly  weakly  umperatura  dapcAdeei.  wiib  Ux 
•>c«{Xloo  of  tbe  TO  pboooo  Aaar  the  X  point  Uoe  Uuct  (c)).  (b)  CalcuUtod  dupenioo  curvea.  FUied  lymbols  ebow  tb«  oooreooraul- 
Ixod  Morflee  Cbar  pbooooi)  and  tbo  open  ctrclee  Indicata  tbc  phoooAi  with  £|  lymmetry,  which  are  renormalLsed  by  ir<U>nuiloe<  wUk 
the  oooductkM  eUarooe.  The  pboooni  with  £4  lyinnKUy  do  not  reoonnellxe  (lAef  Weber  (Ref.  4)). 


Features  of  RU5 


H 

i 


©Provides  the  highest  absolute  accuracy  of  any  routine  elastic  modulus 
measurement  technique. 

©Determines  the  full  ANISOTROPIC  elastic  tensor  in  a  single  measurement. 

©Can  handle  the  smallest  samples  of  any  technique-this  minimizes  sample  prep 
problems  as  well  as  errors  introduced  by  radioactive  heating. 

©LANL  is  the  lead  laboratory  in  the  world  in  the  development  and  use  of  RU5. 
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transducer  system,  and  if  the  sample  is 
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[TR-15]  [TR-16] 


[TR-17] 


RU5  66(56, 6imuItan60U6ly,  the  coilap6e  of  one  oheav  modulue  and  no 
affect  on  the  other  in  a  2mm  6in0le  cry6tal  of  La2.x5r^Cu04  at  the 
6tructural  pha6e  tran6ition  near  223K, 
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[TR-i8] 


;UKJM  siRlKliUn 


Kue  is  a  mature,  fully 
commercialized  technic^ue,  now  in  its 
second  hardware  ^eneration-which 
includes  full  phase  sensitivity  to 
enhance  sharpness 


Resonant  _ 


Appiications  to  fhyii^  Hateriats  Neatunmen 
and  Nondntruct've  Evaluation 


Albert  Miguori 
John  L.  Sarrao 
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Lunof  fockt  Sound  Velocity.  Vp 
ond  cheese*  (Kllofnetert/tecond) 
Lunar  Rock  10017  1.84 

GJetoet  (Norway)  1.83 

Provolone  (Italy)  I.75 

Romano  (Italy)  1,75 

Cheddar  (Vermont)  1.72 

Emmenthal  (Swift)  1.65 

Muentfer  (Wlfconiln)  1.57 

Lunar  Rock  10046  1.25 

(  Science.  168.  1579,  1970  ) 
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The  reeonancee  of  a 
solid,  when  properly 
analyzed  provide  full 
anisotropic  elastic 
information  at  very 
hi^h  accuracy 
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^  .  Following  Hamilton,  we  allow  w,  to  vary  arbit- 

omputstion  of  rosonsness  rarily  in  the  volume  V  and  on  the  surface  S 

(«,  ->  Uj  +  8m,)  and  calculate  the  variation  8L  in 

from  Migliori  et.al.  Physica  B  183,1,1993 


The  "Calvin  and  Hobbes"  model  of  the  vibrations  of  a  rectangular 
parallelepiped 
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ed  crossings  may  be  spun 
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[TR-27] 


How  well  doee  RU5  do? 


f  in  MHz 

mode  f (measured)  f(fit)  %err  wt  k  i 

1  0.340010  0.339508  -0.15  0.00  4  1 

2  0.447200  0.447009  -0.04  1 , 00  6  2 

3  0.490240  0.490304  0.01  1.00  7  2 

4  0.541860  0.541385  -0.09  1.00  3  2 

5  0.557460  0.556670  -0.14  1.00  4  2 

6  0.591740  0.592201  0.08  1.00  5  1 

7  0.608670  0.608409  -0.04  1.00  6  3 

8  0.615590  0.615232  -0.06  1.00  2  2 

9  0.623100  0.622549  -0.09  1.00  8  2 

10  0.642700  0.642863  0.03  1.00  1  2 

11  0.675510  0.676327  0.12  1.00  5  2 

12  0.683560  0.683290  -0.04  1.00  8  3 

13  0.690650  0.690379  -0.04  1.00  2  3 

14  0.746720  .  0.747308  0.081.0053 

15  0.753860  0.754061  0.03  1.00  8  4 

16  0.827430  0.827516  0.011.0073 

17  0.848300  0.849020  0.08  1.00  1  3 

18  0.855320  0.854900  -0.05  1.00  4  3 

19  0.870540  0.870830  0.03  1.00  6  4 

20  0.878390  0.878949  0.06  1.00  3  3 

21  0.882080  0.881934  -0.02  1.00  5  4 

22  0.884870  0.884495  -0.04  1.00  1  4 

23  0.887380  0.887736  0.04  1.00  74 

24  0.891190  0.891837  0.07  1.00  2  4 

rms  errors  0.0653  H 


df/dc,j 

0.00  1 .00 
0.13  0.87 
0.16  0.84 
0.02  0.98 
0.00  1.00 
0.05  0.95 
0.30  0.70 
0.02  0.98 
0.07  0.93 
0.12  0.88 
0.10  0.90 
0.05  0.95 
0.16  0.84 
0.10  0.90 
0.05  0.95 
0.06  0.94 
0.03  0.97 
0.08  0.92 
0.12  0.88 
0.11  0.89 
0.33  0.67 
0.21  0.79 
0.20  0.80 
0.04  0.96 


5120  Steel 

Fine  grained  martensitic 
steel 

free  moduli  are  c11,  c44 
units;  I0'%n/cm* 

mass=0.2875  gm 
0.412  X  0.340  X  0.264  cm 
P=  7.785  gm/cc 
C11=2.7187  +  .19% 
C44=0.8148  ±.02% 


H 
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Ultimate  accuracy  determined  by  geometry-for  a  SijN^  ball 
bearing,  geometry  errors  are  less  that  1  part  in  10®.  So  are  the 
modulus  errors! 
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[TR-35]  [TR-36  Unavailable  at  Time  of  Printing]  [TR-37] 


RUS  can  b©  US6d  at  tsmpsraturss  as  high  Resonances  of  a  5120  steel  RPR  at  38C  and  378C 

as  1800  C  (O.  Anderson  et.  al.).  For  more  measure  using  metal  diffusion  bonded  LiNbOg/Alumina 

moderate  temperatures  (700  C),  both  transducers 


(siui  •qK)epni!fduiv 
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Parallelepiped  sample  held 
lightly  between  transducers 


The  transducers  are  mounted  beneath  the 
0.500"<D  304SS  screens  that  the  sample 
touches.  Only  the  bottom  one  is  visible. 
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RESONANT  ULTRASOUND  SPECTROSCOPY 
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>tain  a  measurement,  but  it  would  be  unnecessarily  costly  to  wait  much 


(rpr.exe)  developed  by  the  authors,  and  the 
workhorse  code  in  their  laboratory,  available  from  the  author  [7.1],  is  a  good 
example  to  use  here  in  expounding  upon  the  nuances  of  obtaining  a 
satisfactory  fit  and.  hence,  moduli.  From  the  spectrum  of  Figure  7.11  we 
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RESONANT  ULTRASOUND  SPECTROSCOPY 
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GInzburg-Landau  description  of  second  order  phase  transitions 
Linear  coupling  between  strain  and  order  parameter  q 
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2.  Compute  the  elastic  response  under  the  constraint. 


‘^44(7')  -  C44  +  producing  a  so-called 

Curie-Weiss  behavior  on  either  side  of  the  transition. 
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(7)  =  3t  quadratic  coupling  between  strain  and  order  parameter 


What  happens  when  coupling  is  turned  on; 


C44(10^^dyn/cm 


RU5  sees,  simultaneously.the  collapse  of  one  shear  modulus  and  no 
effect  on  the  other  In  a  2mm  single  crystal  of  La2-xSr^Cu04  at  the 
structural  phase  transition  near  223K. 
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[TR-64] 


Ocn'AUBOM-  -nUT  rcffW-C>CAy>J-L  ro  e>A.r-nokH0H6ic. 
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[TR-69] 

Electronic  Structure  of  the  Narrow-Gap 
Semiconductor  FeSi  using  RUS 


Elastic  Constants  of  FeSi  Simple  model  for  the  unusual  band-edge 

density  of  electronic  states  in  FeSi 
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1 100  K,  W—'SOO  K,  and  g =4.20  states/cell. 
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JT  a-*  A  oi  raxv/f  iii\  a 


Table  3. 1 -CRYSTAL  STRUCTURE  DATA  FOR  PLUTONIUMt 


Phase 

stability 

Range, 

•c 

Space  Lattice 
and  Space  Group 

Unit  Cell 
Dimensions, 

A 

Atoms  per 
Unit 

Cell 

X-ray 

Density, 

g/cm’ 

Refer¬ 

ence 

a 

Below  ^115 

Simple  monoclinic 
P2,/m 

a  = 

e  21*0: 

6.183  ±  0.001 

b 

C  s 

4.822  ±  0.001 
10.963  *  0.001 

16 

19.86 

6 

13  = 

101.79*  ±  0.01* 

0 

-115  -  -200 

Body-centered 

monocllnlc 

/2/mt 

a  = 

b  = 

c  = 

0  = 

@190*C: 

9.284  *  0.003 
10.463  ±  0.004 
7.859  ±  0.003 
92,13*  ±  0,03* 

34 

17.70 

8 

y 

-200  -  310 

Face-centered 

orthohomblc 

Fddd 

a  = 
b  = 

c  = 

@235*C: 

3.159  ±  0.001 
5,768  ±  0.001 
10.162  ±  0.002 

8 

17.14 

9 

6 

310  -  452 

Face -centered 

@320  *C: 

cubic 

Frriim 

a  = 

4.6371  ±  0.0004 

4 

15.92 

10 

452  -  480 


Body-centered 

tetragonal 

lA/mmm 


@465*0: 

a=  3.34*0.01 
c  =  4.44  ±  0.04 


16.00 


480  -  640 


Body- centered 
cubic 
ImZm 


@490“C: 

a=  3.6361  *  0.0004 


16.51 


’  FTEUln”ger!lctrCrrsr‘‘l6;  "e/aX’  w'h  "iTt’  «•  Zaohariaaen  and 

p>  4.^1  '  j  f  M  ^  '  63)  I  W,  H.  Z^c)i&ri3.SGn  sind  F,  H.  Ellinccr  Acts.  Crvst 

t  A»hn!ih  ”■  AJMETranaartlnna  206:  1256  (1956)  —  ^ 

^nlon  o?cryatXrraDhv”lmn?lr  °/ tabulated  tn  the  International 
oulun  oi  ^..rysiaiiography  International  Tables  fnr  Y-rav  tr_,  , 

1  ^^^^nuTi^UtWs  reta.n^  t^^^TlTag:} 


[TR-77] 


Pu  has  the  largest  shear-wave  anisotropy  of  any  FCC  metal.  This 
can  produce  distortions,  warping  and  a^in^  effects  stronger  than 
expected  for  other  metals— 


ELASTIC  PROPERTIES  OF  FACE-CENTERED-CUBIC 
PLUTONIUM* 

H.  M.  LEDBETTER 

CryofetiKs  Divaion.  Intthitie  ht  Bute  Sund»rdi.  NilioMt  Bureau  of  Sundardi.  Boulder.  CO  JJOMl  U.Sj^ 


R,  L.  MOMErvnr 

Rockwell  Inienitiionit  Rocky  Flan  tttni.  CoMen.  CO  10401.  U.S> 


(RecenW  22  Drcrmbtr  1975) 


\lau  Jcntti): 


I.  tapenmcntal  »uM  »e|odi«  .„d  eUaic  cornunn  for  a  Pu-I  wi  V.  Ga  .juk  cmial  fvefocirie 
_  »  from  <II0>  in  a  |00l)  pbrlei  ^  ‘ 


f  -  15.75  t  O.Ol  t/ci7i‘ 


Loofnudmal  velocity: 
Tnfuvrree  «efocti»c(: 


a,  -  102S4  ±  0.0032  a  lOWt 
a..  -  1.4604  +  a0023  »  lO’cm'i 
t>„  -  0.5674  +  0.0009  a  10*  an/% 


Ebitic  itilTnenca: 


t  -ifC,,  -C,j|-0.47|4  0.032  »  IO'*N  ni' 
C„  -  3.621  I  0.036  a  J0**N/m  ' 

a  -  irr  ^  ^  •  *®'* N/m* 

»  -  KC,,  +  2C,,)  -  2.99J  ±  0.030  a  10“*  N/m‘ 


to 


Elatiic  aiufoiropfos: 


^  "  C'aVC  «.  703  (after  Zetierl 

-  /rX4-l)i  ^  Bueiiem) 


The  shear  anisotropy  in  6- 
stabiiized  Pu  at  room 
temperature  is  7:1,  and  only 
one  measurement  of  it  has 
ever  been  made 
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RUS  has  important  applications  to  nondestructive  testing  for 
quality  control  in  industry  and  government. 


First  launch  of  the  Trident  li  SLBM 


NOT  LIKE  THIS 
(from  The  Economist) 


Resonances  can  be  used  to  determine 
manufacturing  flaws.  Below  is  shown  a 
spectrum  taken  from  a  Si3N4  ball  bearing 
used  in  the  space  shuttle.  The  error  in 
roundness  causes  two  degenerate 
resonances  to  split.  The  splitting 
determines  the  error  in  roundness. 
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ACOUSTIC  RESONATORS  AND  THE  PROPERTIES  OF  GASES 


Michael  Moldover 

National  Institute  of  Standards  and  Technology 
Gaithersburg 
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■~Z  ^  /O  (sf>j/lc  m«f-  r*e>-,  /oyij 


i^dloopci  u 

.P(  C^UiuS 


toltr-f  H»-t:  L  =  "^  jr  “^.J..v) 

^/rchr^tXtiMtrJ  Q  }  EMp  ^ 

Q=^^(g-r^) 

'T 

C*^ri,w^  B-f^t  Jr/—  r^*^' 

'  ^or 

Examples  of  Primary  Thermometers 


Dilute  Gas  ^  5  ^ 

(equation  of  state)  ^  • "  • 


Dilute  Gas  atc*  =  £' 

(speed  of  sound) 


%  =  -V^rci  .  . 


Acoustic  Gas  Thermometry 
Jilute  monatomic  gas 

i.TT  . 


-L  =  -  ^  -  c\n 

T^  273.16  K  c>(rj 
Resonance  Method 


\/ 


*  =  N^k.  =  ltf£l 


xyx  eigenvalue 


.  Mzir”x  f/in]' 


Black  Body  Radiation  2>  .  *>  . 

(power /area).  7  ' 


Johnson  Noise 
(power) 


Dipole  in  Field 
(energy) 


=  41  rAA 


E  =  -fi^B  =  -//B  cfnh 


fe)-(i7. 


Jllo-fi 


Won  0/0-,/i;  5«i  T,m  ^t~J3.l(,^>cad 


±M<v^>=iliX, 


<v">  = 


o/e-fi^^i  fcg 

‘/e-f/'ift  R 


p  -  Mc^  _  /W  fid 

J:  7;^  -^tiBJCooo 


^  o/"  Tut  t-  /f**/^'  ^ 

Volome  erf  r-  If^tf-  ^ 

(D  M 

(B  Tw 
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RAOIAL  RBSONANces  fN  A  st>H£Re 


Ooen  Gas 
and  Vacuum 


Vafva 

Oporator 


W/aHE'5T  PO^SIBLB  Q  AT  Louj 

OBmtTY  Two  VISCOUS 

f/OT  Strive  TO  SMAU.  emoRS 

XU  COUSTRVCTfOU  -  NOU-oecSeUBRATB 

hooes 


Rgure  5.  Cross-section  of  resonator  and  pressure  vessel.  T1 
transducer  assemblies  arc  indicated  by ‘T*.*’  and  the  locations 
the  capsule  thermometers  arc  indicated  by  **PRT/’ 


Simplest  Theory 


Rigid  spherical  shell,  ignore  temperature 


y  =  nth  ipherical  Bessel  f««ctioa 


^  0  .  p  ^  IIL 

da  ^ 


From  Morse  and  Ingard,  Eq.  (9.4.14);  (IP  is  the  velocity  potential) 

...  /C  s 

tJ-ig  =  - - — - - — 

4n  .fHJl(f^d\/ 


For  the  thermal  boundary  layer 

P,  =(1 

c 


(Raytdgh) 

#  of  root;  /  =  0, 1, 2, . .  • 

#  of  Bcssd  fuctiott;  n  =  0, 1, 2, .  •  - 


Upon  integrating  the  bessel  functions,  you  get  the  "right"  answer; 

g  (Y-1)5f  1 

f~  2a  i-n(n  +  1)/z^. 


Note:  Dj 


X 

pCp 


Note:  P,  has  equal  real  and  imaginary  parts;  therefore 
frequency  decrease  -Af  =  g 

One  may  write  down  P’s  for  viscous'boundary  layer,  ducts,  elastic 
response  of  shell,  etc.  (for  radial  and  for  non-radial  modes) 


i 

I 
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Pressure  (MPa) 

FIG.  9.  Measured  resonance  frequencies  minus  calculated  frequencies 
(sealed  by  10*/frcquency)  for  (0^)  modes.  Here,  the  calculation  includes 
the  effect  of  the  thermal  boundary  layer  and  holes  in  the  resonator;  how¬ 
ever,  the  calculation  omits  the  effect  of  shell  motion.  The  linear  dependence 
of  A//  /on  pressure  is  a  result  of  shell  motion.  The  slopes  depend  upon  the 
proximity  of  the  gas  resonances  to  the  shell  breathing  resonance  near  20.2 
kHx-  /  t  A 


of  frequency.  The  points  are  the  average  slopes  of  the  curves  in  Fig.  9.  The 
curve  is  calculated  for  an  isotropic  seamless  shell  using  the  theory  of  elasti¬ 
city  and  the  elastic  constants  tabulated  for  aluminum  (sec  Table  V).  The 
idealized  shell  has  a  breathing  resonance  near  20.2  kHz. 


9/ee/ 

?  hh^ 


Figure  7.  Perturbations  to  the  frequency  and  half-width  of  (0,8) 
mode  as  a  function  of  frequency,  with  argon  in  the  resonator  at 
100  kPa.  The  frequency  was  swept  by  changing  the  temperature 
of  the  resonator,  which  changes  the  speed  of  sound. 
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Table  9.  Sensitivity  of  c©  to  impurities 


Impurity 

M 

To 

1 

ci 

d(ci) 

dx 

(g/mol) 

in  He 

in  Ar 

Hi 

2 

1.4 

0.23 

0.68 

He 

4 

5/3 

0-9 

HiO 

18 

1.32 

-3.93 

0.12 

Nc 

20 

5/3 

*-4,0 

0.5 

Ni 

28 

1.4 

-6.27 

0-03 

Oi 

32 

1.4 

-7.3 

-0.07 

Ar 

40 

5/3 

-9,0 

CO2 

44 

1.4 

-10.3 

-0.37 

Kr 

84 

5/3 

-20.0 

-1.1 

Xc 

131 

5/3 

-31.8 

-2.3 

Hg 

201 

5/3 

-49.0 

-4.0 

Table  10.  Spied  of  sound  ratio  detenninations 

Gas  Comment  njt/’  ^Cgas) 

Vc(Ar-M) 


unprocessed 
purified  26  h 
purified  120  h 
purified  240  h 
purified  240  h 


Vc(Ar-M) 

0.22 
0,27 
0.35 
-191, 5* 
-184.63 
-183.92 
-184.35 
-184.00 


Pressure 

(kPa) 


May  1.  1987 
May  2.  1987 
May  21,  1987 
May  5,  1987 
May  4.  1987 
May  14.  1987 
May  20,  1987 
May  22,  1987 


•  T^c  value  listed  is  the  mean  determined  from  the  (0.2H0  6)  modes  The  rm^  rt  •  r - 

the  mean  for  a  single  ratio  was  I  O  nr.m  r...  ,1.  mooes.  1  he  rms  deviation  from 

other  cases.  ""processed  Ar-40  and  about  0. 1  ppm  for  all 
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Tabic  I.  One-sigma  uncertainties  (in  parts  per  million)  from  var- 
lous  sources  in  the  redetermination  of  R 


J  (Volume)^* 

density  of  mercury  at  20  ’C 
storage  and  handling  of  mercury 
thermal  expansion  of  mercury  (0-20  *C) 
random  error  of  volume  measurements 
corrections  from  weighing  configuration 
to  acoustics  configuration 
mass  of  counterweights 

n  Temperature 

random  error  of  calibrations 
temperature  gradient 


fo.n 

I  0.20 

I  0.67 
/  0.20 


O.IO 

V0.I4 


0.8?  f  0  8 

I  0.4 


III  Af/yo 

Ar-40  standard 

comparison  of  working  gas  to  Ar-40 


O.?/  S  0-7 

I  0.4 


IV 


Zero-pressure  limit  of(/i,/t/f^y 
s.d.  of  Co*  from  70  observations  at  14  pressures 
thermal  boundary  layer  correction  (0.3%  of  ^ 
thermal  conductivity)  T t 

possible  error  in  location  of  transducers 


I 


0.68 

0.30 

0.55 


Square  root  of  the  sum  of  the  squares 


f 

OCL 


A.  *A.,  p-‘ 

/  iOjspH 

t 


5:+c)i£>.o»£5^ 
irs ±o.r  Tke-r, 


^  0^^  4  t . 


[  ^  M  ^  /  M  ye  sg- V  ® 
•  Q— I 


eM:-C 


M  -  —  - 

:x^.u  sgjujz 


J3i  ^ 


Ji  ® 


•CJ»SS: 


i)  *£> 


tCE  I  fipjajctoo 
A  ' 

O  I  <C.  71^  R5  •  Q 

<0  B  •figsstraa 

gST^iil 


siig 

o  *c  T\  :> 
*07C  ^ 
*  •  i 

(O.  ili 
fit  d 
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Figure  20.  Excess  half^widths  of  (0,/i)  resonances  with  argon  in 

the  resonator  scaled  by  lO'^/frcqucncy.  Ag  =  measured  g  minus  Figure  21.  Excess  half-widths  of  (0,n)  resonances  with  helium 

calculated  g.  in  the  resonator  scaled  by  lO^/frequcncy.  Ag=  measured  g 

minus  calculated  g. 


FREQUENCY  (kHz) 


Figure  7.  Perturbations  to  the  frequency  and  half-width  of  (0.8) 
mode  as  a  function  of  frequency,  with  argon  in  the  resonator  at 
100  kPa.  The  frequency  was  swept  by  changing  the  temperature 
of  the  resonator,  which  changes  the  speed  of  sound. 


Figure  8.  Relative  amplitude  of  the  acoustic  pressure  as  a  func- 
tior.  of  frequency  in  the  vicinity  of  the  (0,2)  and  (0,7)  modes. 
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■.««  0,.  d»cu>,  Z  ' 

.i.~db>d  .rij  toodoTta  Xl.tTJr  *" 

resonances  at  4777.63  and  4779  38  Hr  f  *Sht  parameters  spec 
«n«r  background  ..mb.  £  M  -  '  * 

fitted  function.  TTic  fitted  u  •  Measured  voltages  min 

ctersspccify  resonances  at  4777  693  ^777^3  ^1^477^  P®''®' 

widths  of 0.550.  0.546  and  0  555  i  S  ’  ^a 

cify  i?and  Cin  Eq.  (75^  parameters  sp 
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0  10,000  20.000 


Frequency  (Hz) 

FIG.  16.  Measured  zcro<pressure  frequencies  minus  calculated  frequencies 
for  the  (1.S)  modes.  The  solid  symbols  represent  the  zero  pressure  inter¬ 
cepts  of  straight  lin«  fitted  to  data  such  as  those  displayed  in  Fig.  10.  The 
curves  arc  obtained  from  Eqs.  (68)  and  (69)  with  the  parameters 
€o  =  3.5  X  10”“*  and  e,  3. 1 X  lO”**.  The  open  symbols  arc  averages  of  the 
three  zero-pressure  frequencies  for  each  ( l.s)  set  of  modes.  These  averages 
are  also  plotted  in  Fig.  15. 


Frequency  (Hz) 

FIG.  15.  Measured  zero-pressure  frequencies  minus  calculated  frequencies 
for  the  (0.r)  (solid  syml^Is)  and  (1.5)  (open  symbols)  modes.  The  sym¬ 
bols  represent  the  zero  pressure  intercepts  of  straight  lines  fitted  to  data 
such  as  those  in  Figs.  9  and  10,  The  intercept  for  the  (0,9)  mode  at  21.4  kHz 
is  254  parts  in  10*  above  the  predicted  value.  The  (0,9)  mode  is  close  to  the 
resonances  in  the  vent  hole  and  coupling  duct. 


Pressure  (MPa) 

FIG.  10.  Measured  resonance  frequencies  minus  calculated  frequencies 
(scaled  by  10*/frcqucncy)  for  the  three  components  of  the  (1.2)  and  (1,6) 
resonances.  Here,  the  calculation  includes  the  effect  of  the  viscous  and  ther¬ 
mal  boundary  layers;  however,  neither  the  effects  of  shell  motion  nor  the 
effects  of  the  holes  arc  included. 


*'«Ponsc  of  the  shell  to  exciution  with  symmetry  of 
(0,^)  as  a  function  of  frequency.  The  points  arc  the  average  slopes  of 
eulves  such  as  those  shown  in  Fig.  10.  The  curve  is  calculated  for  an  isotrop¬ 
ic.  seamless  shell  using  the  theory  of  elasticity  and  the  clastic  constants  tabu  - 
latcd  for  aluminum.  The  idealized  shell  has  resonances  at  0  and  24.5 
kHz. 


J 

1 
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AMPLITUDE  DEVIATION 


Measurement  of  the  ratio  of  the  speed  of  sound  to  the  speed  of  light 

James  B.  Mch! 

Physics  Department,  University  of  Delaware,  Newark,  Delaware  19716 


Michael  R.  Moldovcr 

Thermophysics  Diuision,  National  Bureau  of  Standards.  Gaithersburg.  Maryland  20899 
(Received  2  June  1986) 

Mewurements  of  the  resonance  frequencies  of  the  acoustic  modes  and  of  the  microwave  modes  of 
single  cavity  can  determine  u/c.  the  ratio  of  the  speed  of  sound  of  a  gas  to  the  speed  of  ligh/ 

unnij^“,T'"‘*  a  monatomic  gas  would  determine  the  thermodynamic  temperature  T  with 
nprec^ented  accuracy.  By  judicious  choices  of  cavity  geometry  and  resonance  modes,  u/c  can  be 
th^r^  Part-per-milhon  accuracy  using  cavities  whose  geometry  is  known  only  to  parts  per 
thousand.  These  techniques  can  also  be  applied  to  measurements  of  the  universal  g«  consunt^ 

"he m“gl  "  average^tomicr  of 


<473,5  1474.0  1474.5 


MODE 


6  V(V-V(V 
^  V(n) 


yALF  it 

-U//0TAj 


TM/(  If  is.  I*  0.7  fiHeJ  s 

IflS.S-d  f.O  caU  g 

ifn.oto.f  t.cic  I 

TMii  i^is.Sr*  o.e  h'TieJ  <j 

iffO.I  ±  <Xe>  Calc  9 

Hie-2±l.f  cole  I 


FREQUENCY <MHz)  MERCOi.'/  /4  /  g  6  ±  /  5" 
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r/K 


» 


106x 


theriwodywaiwic  properties 

1.  Ideal-gas  heat-capacity:  C^(T] 

2.  Virlal  equation  of  state: 

pV=RT{l  *...) 

Note:  Resonance  techniques  are  not  recommended 
for  liquids  because  oscillations  of  container  cannot 
cannot  be  separated  from  oscillations  of  fluid 
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Spe«d  of  sound  dau  for  CHr,<CF,. 
Psntanuoroethsno  -  refrigerant  R12S 

Temperature  range  0.71  c  TiT^t,  1.18 


Oevfatlotts  of  dau  from  o</*«7)  surface. 
RMS  deviations  are  ZxtO^Ki/. 

Impurities  cause  larger  systemacic  errors. 
The  surface  uses  3  parameters  for  C,*(71. 
fideal-oas  heat<apadt^.  3  parameters 
for  each  of  the  virial  coeffictents  8f  7)  and 
C(7)  aiMl  2  parameters  lor  D{T\. 


Ideal  gas  heat  capacity  of  Rl  25. 
deduced  from  the  speed  of  sound  data 
through  the  relation: 

-  Mu/ntAu,^^  RT\. 

/f  is  itw  universal  gas  constant;  M  is  the 
•nolecula^  weight  and  u/  is  the  jero- 
pressure  fimit  of  i/V.  7) 


Oeviatioos  of  C^f  HW  from  a  polynomial 
o  present  dau 
■  independent  spectroscopic 
information  (estimatedi 


Calibration  with  argon;  we  measure 

method  requires  stable  resonator,  frequency  standard, 
and  thermometer;  however,  many  errors  in  calibration 
can  be  tolerated. 


tolerant  of  temperature  gradients  in  bath 

example:  propane,  temperature  range  210  K  -  460  K 
^propane  changes  ^05% 

(^afxoo^^propaae)  ChangGS  16% 


C^(T) 


1 


Mu^  -  RT 


3  M 

\  ~  -1—1 

5  M 

p 


'/Y 

J.) 


to  obtain  Cp®(7)  to  0.1  %  requires  T  to  0.23  K 

MJM^  to  0.0001 

intolerant  of  ImpuritieSr  vexing  problem  with  mixtures 

inconsistencies  among  modes  <  0.0025% 


compound 


lemperature  maximum  number  of 

range  pressure  isodierms 

K  kPa 


CaadidaCc  refrigeranfs  for  vapor  compression  cycles 


CF,-CHFj 

RlMa 

233  -340 

600 

10 

CFOj-CH, 

RMlb 

260  -  315 

.70 

5 

CHa,-CF, 

R123 

260-  335 

80 

6 

CHFa-CF,-CI 

R123a 

265  -  300 

50 

2 

CHF,-0-CHFj 

EI34 

255  -  327 

170 

6 

CHF,-0-CHF, 

E134 

255  -  374 

90 

8 

CF,-0-CH,-CF, 

E245 

278  -  384 

50 

5 

CF,-HF-CHF, 

R236ca 

267  -  380 

600 

8 

CHFa-CF, 

R124 

250  -  400 

900 

17 

CHF,-CF, 

R125 

240-400 

1.000 

9 

CHF,-CH, 

R152a 

240  -  400 

l.OOO 

9 

CF,-CH, 

R143a 

240-400 

1. 000 

9 

CF,-CH,-CF, 

R236fa 

276  -  400 

1,000 

7 

CHF,-CF,-CH,F 

R245ca 

311  -400 

900 

5 

CF,-0-CF,H 

EI25 

260  -  400 

1.000 

13 

1  composition 

RI34a/R32/R125 

260-400 

1,000 

12 

CF,-CF,-CF,-CH,F 

R338mccq 

300  -  400 

400 

6 

*  Thermoacoustic  Refrigeration 

5  compositions 

HenCe 

210-400 

1.500 

42 

Semiconductor  Processing 

SF. 

230  -460  1.500 

16 

CF, 

300  -  475  1.500 

9 

C,F* 

210  -  475  1.500 

14 

Cl, 

HBr 

BCl, 

WF. 

planned 

! 

correlation  of  u(p,T) 


<  0.002%  r.m.s 


excess  half-widths:  t^glf  -  0.004%  -  0.02%, 

depending  upon  mode 
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TRAWSPORT  PROPERTIES 


1.  Greenspan  acoustic  viscometer  for  z;; 


/Ve//t 


(mkT)' 


hard  sphere  of  diameter  o 


as 

1  10 


2.  Prandtl  number  machine 

pr  ~  ^  -  viscous  diffusivity 

A  thermal  diffusivity 

=  —  for  hard  sphere 
3 

3.  If  time  permits,  electromagnetic  equivalent  of 
Greenspan  viscometer:  reentrant  resonator 


<j 

o  0 


10 


100 

TEMPERATURE/  K 


I  (mm) 


Loss  mechanisms  spatially  separated  in  standing  waves. 
High  thermal  losses  High  thermal  losses 


t 

High  viscous  losses 


Greenspan  viscometer:  A  length 


Double  Helmholtz  Resonator 


Wo 


,  ^r_L 

L  iv,  V4.J 


UJ 


•  i  )  (i 

ItO-cXr-O  ^ 


'ey 

Acoustic  pressure  contours  rfear  duct  end 
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incidence  (solid  lines)  tor  each  of  the  three  models  discussed. 
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Figure  3.4.  Acoustic  phase  velocity  measured  with  the  probe  microphone 
In  a  sand  with  a  flow  resistivity  of  85  rayls/cm- 
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PHYSICAL  PROPERTIES  OF  POROUS  MATERIALS 
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SIMPLE  POROUS  MEDIA  MODEL 


T|  =  coefficient  of  viscosity 
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Bessel’s  equation  and  solution  for  velocity 


f)u.  in  u  cylindrical  tube  having  a  circular  cross-section,  as  a  function  of  s. 
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ACOUSTIC  INVERSION  FOR  PORE  PROPERTIES 
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COUPLING  OF  SOUND  INTO  PORO-ELASTIC  SOIL 
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velocity  and  attenuation  increase  with  increasing  frequency 
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CONSIDER  SOLUTIONS  FOR  WAVE  EQUATIONS 
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Loudspeaker  source 


University  of  Mississippi 


'  Loudspeaker  source. 

•  Travel  time  curve  using  a  1  kHz  tone  burst. 

•  Porous  medium  is  composed  of  commericially  available  sandblasting  sand 

•  The  signal  on  probe  microphone,  in-situ  microphones  and  geophones  are 
dominated  by  the  type  II  P-wave  at  depths  less  than  about  1 5  cm. 

•  The  type  I  P-wave  dominates  at  depths  greater  than  20  cm. 


Shaker  source 


2D  30 
Distance  (cm) 


f  =  750  Hz 
r  s  1  kHz 


The 

University  of  Mississippi 


Mechanical  shaker  source. 

Travel  time  curves  using  3  different  tone  bursts. 

No  measurements  with  the  probe  microphone. 

The  in-situ  microphones  and  geophones  are  responding  to  the  same  wave. 
The  single  slope  indicates  the  presence  of  only  one  wave  type,  i.e.  Type  I 
P-wave. 

'  No  dispersion  in  velocity  over  this  frequency  range. 

-  The  measure  velocity  is  the  same  as  measured  at  depths  greater  than  20  cm 
using  the  loudspeaker  source. 
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Figure  A4.  The  influence  of  Increasing  the  wave-speeds  In  the  substrata 
factor  of  10^  on  the  predicted  normal  surface  impedance  o 
poro-elastlc  layer  at  normal  Incidence  (solid  line)  an 
(dashed  line)  Incidence. 
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